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ABSTRACT 


A  wind-tunnel  study  to  find  the  lift  and  drag 
characteristics  of  a  low-aspect-ratio  wing/body  configuration 
from  an  angle  of  attack  (AOA)  of  -8  to  50  degrees  was 
conducted.  A  further  study  to  find  the  comparative  lift 
enhancement  using  the  same  wing/body  with  a  close-coupled 
canard  for  wing/body  angles  of  attack  of  10,  22,  34,  40,  and 
48  degrees  and  canard  deflection  angles  from  -25  to  25  degrees 
was  carried  out.  It  was  found  that  a  properly-located  canard 
enhanced  the  lift  at  all  tested  angles  of  attack,  compared  to 
the  baseline  wing/body  configuration  results.  The  lift 
enhancement  was  maximized  in  the  post-stall  regimes,  reaching 
values  up  to  34%.  A  small  improvement  in  lift-to-drag  ratio 
was  noted  at  all  tested  angles  above  10  degrees  angle  of 
attack. 
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I .  BACKGROUND 


A.  AIRCRAFT  AGILITY 

In  future  air-to-air  combat  scenarios  the  aircraft  that  is 
the  most  maneuverable  while  at  the  same  time  controllable  will 
in  all  likelihood  win  the  engagement.  Winning  an  air-to-air 
engagement  against  another  agile  fighter  using  all-aspect 
missiles,  such  as  the  AIM-9L,  requires  the  ability  to  out 
maneuver  the  opposition  to  point  and  shoot  first.  [Ref.  1] 

In  close-in  combat  with  another  agile  fighter  the  aircraft 
that  "wins"  will  have  three  advantages  over  its  opponent. 
First,  the  aircraft  will  need  the  ability  to  "outpoint"  its 
opponent.  The  aircraft  that  points  at  its  opponent  first  has 
the  first  opportunity  to  launch  a  "fire-and-forget"  weapon. 
Second,  the  aircraft  must  have  the  ability  to  keep  the  nose  of 
the  aircraft  pointed  at  the  opposition  for  a  longer  time.  The 
aircraft  that  can  point  longer  has  the  capability  to  maneuver 
at  higher  turn  rates  for  longer  periods  than  the  opposition. 
Maintaining  nose-point  longer  than  the  opposition  enables  the 
aircraft  to  defend  itself  against  aircraft  other  than  the 
target  or  to  engage  multiple  targets.  Third,  the  aircraft 
must  have  a  greater  straight-ahead  acceleration  than  the 
opposition.  Greater  acceleration  gives  the  ability  to  escape 


1 


the  battle  or  to  reengage  the  opposition  to  the  advantage  of 
the  agile  aircraft.  [Ref.  1] 

To  outpoint,  turn  faster,  and  accelerate  faster  than  an 
opponent  the  aircraft  must  be  supermaneuverable . 
Supermaneuverability  is  a  blend  of  post-stall  maneuvers  and 
the  use  of  sideslipping  or  direct  force  control.  Post-stall 
manueuverablity  (PST)  is  the  capability  to  tactically  maneuver 
the  aircraft  in  a  controlled  fashion  beyond  the  maximum  lift 
angle  of  attack.  Direct  force-control  maneuver  (DFM)  is  the 
ability  to  change  the  aircraft's  yaw  and  pitch  independent  of 
flight  path  or  to  maneuver  the  aircraft  at  constant  fuselage 
attitude.  A  PST  maneuver  sacrifices  energy  for  a  positional 
advantage.  A  typical  PST  maneuver  involves  rapidly  pitching 
the  aircraft  to  a  high  angle  of  attack  and  maintaining  this 
condition  for  2  to  3  seconds.  The  aircraft  that  maneuvers  in 
the  PST  regime  will  complete  a  tactical  maneuver  in  less  time 
and  space  than  an  aircraft  not  executing  a  PST  maneuver. 
Unfortunately  there  are  a  number  of  limitations  which  prevent 
a  conventional  aircraft  from  performing  a  PST  maneuver.  [Ref. 
2] 

At  high  angles  of  attack,  the  flow  tends  to  separate  from 
the  wing  surface.  The  center  of  pressure  and  neutral  point 
move  aft  with  respect  to  the  center  of  gravity  as  the  airflow 
separates  from  the  wing  surface  at  the  maximum  lift  angle  of 
attack.  The  control  surfaces  are  also  losing  their 
effectiveness  at  the  same  time.  The  result  is  that  the 
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aircraft  encounters  severe  trim  and  stability  problems  when  it 
has  the  greatest  control  power  requirements  and  is  the  least 
controllable.  At  angles  of  attack  greater  than  30  degrees,  an 
additional  reaction  control  system  is  necessary  for  pitch  and 
yaw.  The  deflection  of  the  engine  thrust  is  one  of  the 
proposed  methods  for  a  reaction  control  system.  In  general, 
at  angles  of  attack  greater  than  50  degrees,  engine  thrust 
vice  aerodynamic  lift  provides  the  necessary  lift  and  control. 
[Ref.  2] 

A  DFM  is  used  to  aim  the  fuselage  at  the  opponent 
independent  of  flight  path  for  more  precise  firing  solutions. 
Studies  have  shown  that  for  a  rear  hemisphere  attack,  an 
elevation  aiming  ability  of  at  least  6  degrees  and  an  azimuth 
aiming  ability  of  about  4  degrees  would  be  advantageous .  With 
conventional  aircraft,  elevation  aiming  can  be  achieved  to  a 
limited  extent  through  the  use  of  wing  flaps  in  conjunction 
with  elevator  deflection.  Azimuth  aiming  for  a  conventional 
aircraft  would  involve  the  use  of  special  control  surfaces  and 
an  associated  flight  control  system.  [Ref.  2] 

An  example  of  a  proposed  supermaneuverable  aircraft  is  the 
X-31.  The  X-31  aircraft  uses  a  long-coupled  canard  for 
controllability  with  thrust  vectoring  to  be  supermaneuverable 
in  the  post-stall  flight  regime.  Figure  1  shows  the 
difference  between  a  close-  and  long-coupled  canard.  More 
will  be  noted  on  the  difference  directly.  When  the  angle  of 
attack  is  increased  to  the  point  where  the  aircraft  begins  to 
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stall,  computer-controlled  thrust  vectoring  comes  into  play. 
The  thrust  vectoring  is  provided  by  three  curved  paddles  that 
can  move  into  the  exhaust  stream  to  deflect  the  flow.  The 
ability  to  successfully  use  thrust  vectoring  assumes  adequate 
control  to  prevent  spin  departure.  [Ref.  3] 


Figure  1.  Types  of  Canards  [Ref.  6] 


Herbst  states,  "Unfortunately  there  is  a  progressive 
effect  of  thrust  versus  weight  on  weight  and  cost  and  even  an 
asymptotic  limit."  [Ref.  21  A  prohibitive  cost  limit  in  the 
race  to  achieve  greater  thrust-to-weight  would  imply  aircraft 
such  as  the  X-31  that  extensively  use  thrust  vectoring  may  not 
be  built  in  sufficient  numbers  to  counter  the  opposition's 
agile  aircraft.  Accordingly,  other  methods  of  obtaining 
maneuverability  and  control  in  the  post-stall  regime  warrant 
increased  investigation.  [Ref.  2) 
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One  such  area  of  increased  interest  is  dynamic  stall. 
Dynamic  stall  is  a  complex  event  where  there  is  a  "...dynamic 
delay  of  stall  on  a  wing's  airfoil  experiencing  unsteady 
motion,  to  angles  significantly  beyond  the  static  stall 
angle."  [Ref.  4]  The  delay  of  stall  is  usually  followed  by 
large  changes  in  lift  and  pitching  moment.  Dynamic  stall  was 
first  observed  on  helicopter  blades.  It  was  found  that 
increased  lift  could  be  obtained  by  rapidly  pitching  the 
airfoils.  The  increased  lift  was  due  in  part  to  the  vortex 
formed  on  the  airfoil  during  the  unsteady  motion.  It  was 
found  that  the  increased  lift  due  to  dynamic  stall  could  also 
be  utilized  by  fighter  aircraft  when  the  aircraft  was  rapidly 
pitched.  This  increase  in  lift  due  to  dynamic  stall,  as 
noted,  is  an  unsteady  time-dependent  phenomenon.  At  the 
present  time  predictable  control  of  this  unsteady,  separated 
flow  for  increased  lift  has  not  been  achieved.  Accordingly, 
other  potential  methods  for  lift  enhancement  such  as  the  use 
of  a  close-coupled  canard/wing  are  receiving  increased 
attention.  [Ref.  4] 

B.  CANARD /WING  INTERACTION 

Increased  agility  through  the  use  of  a  close-coupled 
canard  configuration  for  enhanced  lift  has  been  the  subject  of 
growing  scientific  interest  and  practical  aeronautical 
application.  In  the  1960's  Behrbohm  established  that  a  close- 
coupled  canard  with  the  canard  and  wing  of  aspect  ratios 
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between  1  and  3  had  an  increased  coefficient  of  lift  over 
their  non-canard  counterparts.  The  increased  maximum  Cl  in 
turn  contributed  to  an  increase  in  the  angle  of  attack  that 
could  be  achieved.  The  increase  in  maximum  Cl  was  due  to 
constructive  interference  between  the  vortex  systems  of  the 
wing  and  canard.  It  is  thought  that  constructive  interference 
occurs  when  the  downwash  from  the  canard  suppresses  the  flow 
separation  on  the  wing.  The  formation  of  wing  leading-edge 
vortices  are  delayed  until  induced  downwash  of  the  canard 
supports  flow  separation.  The  longer  the  flow  separation  is 
delayed,  the  greater  the  lift  enhancement.  [Ref.  5] 

The  translation  of  theory  into  reality  was  realized  with 
the  SAAB  AJ-37  Viggen  aircraft.  The  Viggen  aircraft  uses  a 
close-coupled  canard  that  was  able  to  generate  a  65  percent 
greater  maximum  Cl  at  approach  than  a  pure  delta  wing.  The 
Viggen  aircraft  successfully  took  advantage  of  favorable 
vortex  wing-canard  interactions  to  achieve  greater  lift  which 
allowed  the  aircraft  to  lower  its  landing  and  takeoff  speed 
for  a  STOL  capability.  [Ref.  6]  Figure  2  illustrates  the 
canard/ wing  vortex  interaction  of  the  SAAB  Viggen  [Ref.  7] 
The  use  of  a  close-coupled  canard  on  the  SAAB  Viggen  also 
gave  the  aircraft  much  greater  trim  control  compared  to  a 
tailless  delta-winged  aircraft  such  as  the  Mirage  III.  The 
lack  of  a  horizontal  tail  on  aircraft  such  as  the  Mirage  III 
requires  the  use  of  elevons,  which  are  deflected  upward  to 
create  a  down  thrust  to  rotate  the  aircraft  for  landing  and 
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takeoff.  The  elevens  have  a  much  shorter  moment  arm  than  a 
horizontal  tail,  which  in  turn  requires  that  the  elevens  be 
large  to  be  effective.  The  elevens  exact  a  weight  penalty  of 
as  much  as  two  tons,  but  more  importantly  decrease  the 
effective  lift  of  the  aircraft,  just  the  opposite  of  what  is 
needed  for  enhanced  takeoff  and  landing  ability.  The  SAAB 

% 

Viggen  uses  a  large  canard  with  trailing-edge  flaps.  During 

*  takeoffs  and  landings  the  Viggen  has  lift  from  the  main  wing, 

lift  from  the  canard,  plus  lift  from  the  downward-depressed 
main  wing  elevens.  The  result  is  that  the  Viggen  has  a  much 
shorter  takeoff  and  landing  distance  than  the  Mirage  III. 
[Ref.  7] 


Figure  2.  Canard/Wing  Vortex  Interaction  [Ref.  7] 


Stoll  and  Koenig  demonstrated  that  the  maximum  lift  of  a 
close-coupled  canard  model  was  34  percent  greater  than  a  non- 


7 


canard  version  of  the  same  model.  Furthermore,  the  increase 
in  lift  could  not  be  solely  attributed  to  an  increase  in  wing 
reference  area.  The  wing  reference  area  of  the  canard-wing 
model  was  only  15  percent  greater  than  that  of  the  wing-alone 
model.  [Ref.  8] 

Er-El  and  Seginer  found  that  a  close-coupled  canard  placed 
upstream  and  above  a  60-degree  swept  wing  delayed  the  onset  of 
wing-leading-edge  vortex  breakdown  for  an  angle-of-attack 
range  from  14  to  24  degrees.  But  significantly,  the  use  of  a 
close-coupled  canard/wing  did  not  always  result  in  improved 
aerodynamic  characteristics.  Improved  lift  was  dependent  upon 
proper  longitudinal  positioning  of  the  canard.  Er-El  and 
Seginer  did  not  rigorously  establish  what  this  longitudineil 
position  should  be,  but  that  was  not  the  objective  of  their 
study.  [Ref.  9] 

Calarese  tested  a  model  with  the  canard  placed  above  the 
wing  and  coplanar  to  the  wing.  At  all  angles  of  attack,  the 
canard  mounted  above  the  wing  created  a  more  favorable 
interaction  between  the  leading-edge  vortices  than  the 
coplanar  canard.  The  use  of  a  canard  placed  above  the  wing 
caused  a  noticeable  improvement  in  the  lift-to-drag  ratio  for 
the  tested  angles  of  attack  of  10,  16,  and  19  degrees.  The 
increase  in  the  lift-to-drag  ratio  was  12  percent  greater  than 
the  lift-to-drag  ratio  for  the  coplanar  model.  It  was 
surmised  that  the  more  favorable  lift-to-drag  ratio  of  the 
model  with  the  canard  positioned  above  the  wing  was  due  in 
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part  to  the  vortex  and  wakeflow  of  the  canard  missing  the  wing 
and  therefore  causing  less  interference.  A  more  likely 
explanation  is  that  downwash  suppressed  the  flow  separation  of 
the  wing  and  thus  improved  lift.  [Ref.  10] 

O'Leary  and  Weir  demonstrated  that  the  maximum  Cj,  using  a 
close-coupled  canard  was  20  percent  greater  and  at  an  angle  of 
attack  5  degrees  greater  than  for  a  non-canard  model .  At 
angles  of  attack  from  0  to  36  degrees,  the  canard  had  little 
effect  on  the  slope  of  the  lift  versus  angle-of-attack  curve 
for  canard  deflections  of  0,-10, -20  degrees.  A  slightly 
higher  value  for  maximum  Cj,  was  achieved  with  a  canard 
deflection  of  -20  degrees  than  with  canard  deflections  of  0 
and  -10  degrees.  A  positive  canard  deflection  of  10  degrees 
reduced  the  slope  of  the  lift  curve.  At  a  canard  deflection 
of  -40  degrees,  the  maximum  Cj,  was  reduced  by  approximately  20 
percent.  [Ref.  11] 

The  research  by  David  W.  Lacey  at  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center  (now  the  David  W. 
Taylor  Research  Center)  in  the  1970' s  in  the  area  of  close- 
coupled  canard/wing  interaction  was  quite  comprehensive. 
Lacey  found  that  the  increase  in  the  maximum  Cj,  was  mainly  a 
function  of  the  ratio  of  the  canard  area  to  the  wing  area 
(S^/SJ  and  the  canard  placement  in  both  the  longitudinal  and 
vertical  directions.  The  longitudinal  position  was  measured 
from  the  0.25  mean  aerodynamic  chord  of  the  wing  to  the  0.40 
exposed  root  chord  of  the  canard.  For  lift  enhancement,  it 
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was  found  that  the  ratio  of  the  longitudinal  canard  position 
to  the  mean  aerodynamic  chord  of  the  wing  with  respect  to  the 
centerline  of  the  fuselage,  x/Caac»  should  be  between  1.0  to 
1.25.  If  the  canard  was  mounted  too  far  forward,  the 
canard/wing  combination  would  actually  generate  less  lift  than 
what  the  wing  and  canard  could  generate  separately.  (  It 
should  be  noted  that  the  long-coupled  canard  on  the  X-31  is  a 
control,  and  not  a  lifting,  device.)  It  was  determined  that 
the  canard  should  be  positioned  so  that  the  ratio  of  the 
vertical  distance  from  the  plane  of  the  wing  to  the  mean 
aerodynamic  chord,  2/c„ae,  should  be  equal  to  0.2.  The  trailing 
edge  of  the  canard  and  the  leading  edge  of  the  wing  should 
never  overlap,  for  a  loss  of  lift  results.  The  combination  of 
a  proper  selection  of  vertical  and  longitudinal  position 
maximized  Ci,  and  L/D.  Increasing  the  size  of  the  canard 
increased  lift  in  a  fairly  linear  fashion  up  to  S^/S^  =  0.25. 
Beyond  0.25  a  sharp  drop  off  in  canard/wing  lift  effectiveness 
resulted.  [Ref.  6] 

Lacey  tested  canards  with  a  leading-edge  sweep  of  25,  45 
and  60  degrees.  Maximum  lift  was  developed  with  the  60-degree 
swept  delta  canard.  Maximum  L/D  was  developed  with  the  25- 
degree  high-aspect-ratio  canard.  A  tradeoff  study  was 
conducted  and  found  that  for  maximum  lift  and  L/D  the  45- 
degree  canard  should  be  chosen,  closely  followed  by  the  60- 
degree  canard.  If  any  of  the  planforms  were  properly  located 
they  would  enhance  lift.  [Ref.  6] 
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Lacey  noted  the  SAAB  TN  60  report  postulated  that  a 
highly-swept  canard  delta  wing  should  be  chosen  to  maximize 
the  synergism  between  the  canard  and  wing  vortices  for 
increased  lift.  Lacey  found  that  a  45-degree  canard  with  a 
25-degree  swept  wing  also  greatly  enhanced  lift,  but  in  this 
case  neither  the  wing  nor  the  canard  generated  the  strong 
leading-edge  vortices  "required”  by  the  SAAB  report.  Lacey 
postulated  that  the  canard  downwash  delayed  the  leading-edge 
stall  of  the  wing  in  a  manner  similar  to  that  of  a  leading- 
edge  slot.  The  canard  could  be  thought  of  as  a  large  low-drag 
boundary-layer  device.  [Ref.  6] 

It  was  shown  that,  with  the  canard  in  a  high  position 
above  the  wing,  the  overall  lift  was  actually  less  than  the 
wing/body  alone  configuration  for  angles  of  attack  below  18 
degrees.  The  interference  that  took  place  between  the 
flowfields  of  the  canard  and  wing  that  caused  this  lift 
degradation  was  not  very  well  understood.  It  was  thought  that 
destructive  interference  occurred  when  the  wash  of  the  canard 
impinged  upon  the  wing  in  an  upward  direction,  which  then 
promoted  flow  separation  on  the  wing  and  thereby  caused  a  loss 
of  lift.  Beyond  an  18  degree  angle  of  attack,  the  wash  of  the 
canard  impinged  upon  the  wing  in  a  downward  direction,  and 
thereby  delayed  the  onset  of  flow  separation  on  the  wing  and 
correspondingly  increased  the  lift.  [Ref.  6] 

Deflecting  the  canard  in  5-degree  increments  from  -10  to 
+10  degrees  with  the  model  at  a  fixed  angle  of  attack  resulted 
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in  a  change  of  Cl  of  only  0.08.  Positive  canard  deflections 
for  a  main  body  angle  of  attack  of  5  degrees  were  found  to 
increase  drag  and  correspondingly  decrease  the  maximum  L/D 
markedly.  Slightly  negative  canard  deflections  increased  the 
maximum  L/D  because  a  0-degree  local  angle  of  attack  for  the 
canard  minimized  induced  drag.  (Ref.  6] 

C.  OSCILIATIMG  CANARD 

The  next  area  of  increased  interest  for  lift  enhancement 
involves  interactions  between  an  oscillating  close-coupled 
canard  and  the  main  wing  flowfield.  Thus  far  few  studies  have 
been  completed  in  this  area.  Buyer  and  Luttges  investigated 
the  flowfield  interaction  between  the  unsteady  wake  of  an 
oscillating  canard  upstream  of  a  static  wing.  An  NACA  0015 
airfoil  was  used  for  both  the  canard  and  the  main  wing.  The 
main  wing  was  mounted  coplanar  to  the  canard  and  0.5  chord 
lengths  downstream.  Main  wing  angles  of  attack  of  10  and  20 
degrees  were  used.  The  mean  canard  deflection  angle  was  15 
degrees  and  the  oscillation  amplitude  was  +/-  10  degrees.  The 
canard  was  oscillated  about  the  quarter  chord  with  periods  of 
156  and  105  msec.  Buyer  and  Luttges  found  that  the  dynamic 
stall  vortex  from  the  oscillating  canard  energized  the 
boundary  layer  of  the  main  wing  which  resulted  in  flow 
reattachment  at  angles  of  attack  far  exceeding  static  stall 
angles  for  the  main  wing.  But  the  amount  of  enhanced  lift  was 


12 


not  well  quantified.  No  comparison  was  made  to  a  static 
canard  and  main  wing  case.  [Ref.  12] 

Other  studies  have  considered  the  flowfield  due  to  an 
oscillating  canard  for  an  X-29  aircraft  model,  but  no  lift- 
enhancement  results  have  been  presented.  [Ref.  13  and  Ref.  14] 

D.  STATEMENT  OF  PURPOSE 

From  the  previous  work  of  others  it  is  fairly  well 
understood  that  locating  the  canard  in  an  optimum  position 
vertically  and  longitudinally  will  generate  vortices  that  will 
constructively  interfere  for  maximum  lift.  Studies  have  not 
been  conducted  for  main  wing/body  angles  of  attack  greater 
than  36  degrees,  nor  have  static  canard  deflection  angle 
increments  of  less  than  5  degrees  been  used. 

It  is  known  that  oscillating  a  coplanar  canard  separated 
from  the  main  wing  by  0.5  chords  will  reattach  the  flow  over 
the  main  wing  at  angles  far  exceeding  the  normal  stall  angle 
of  attack  for  the  main  wing  alone.  The  canard,  however,  was 
not  optimally  located  vertically  and  horizontally  to  maximize 
lift  enhancement,  and  only  two  main  wing  body  angles  of  attack 
of  10  and  20  degrees  were  used.  Additionally,  the  lift 
enhancement  with  the  use  of  an  oscillating  canard  was  not 
compared  to  the  lift  enhancement  with  a  static  canard.  The 
objectives  of  this  investigation  were  to  conduct  : 
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(1)  A  baseline  study  to  find  the  lift  and  drag 
characteristics  of  a  low-aspect-ratio  wing/body  model  from  low 


angles  of  attack  to  angles  of  attack  beyond  36  degrees. 

(2)  A  further  study  to  find  the  comparative  lift 
enhancement  using  the  same  wing/body  and  a  close-coupled 
canard  for  wing/body  angles  of  attack  of  10,  22,  34,  40,  and 
48  degrees  and  canard  deflection  angles  -25  to  25  degrees. 

Upon  completion  of  the  above  investigation  a  follow-on 
study  will  be  conducted  to  examine  the  lift  enhancement  of  an 
optimally-located  oscillating  canard  for  comparison  to  the 
results  of  the  static  canard/wing  configuration. 


14 


II.  EXPERIMENT  AND  PROCEDURES 


A.  OVERVIEW 

A  close-coupled  canard  model  was  designed  and  constructed 
at  the  Naval  Postgraduate  School.  A  calibration  rig  was  then 
designed  and  constructed  to  calibrate  an  existing  external 
strain-gage  balance  in  the  NPS  low-speed  wind  tunnel.  Data 
were  acquired  from  the  signal  conditioning  assembly  through  a 
digital  multimeter,  multiplexer  and  amplifier,  and  stored  on 
a  floppy  disk.  The  results  were  reduced  to  lift  and  drag 
coefficients . 

In  a  baseline  run  of  the  model  without  the  canard,  the 
angle  of  attack  of  the  wing  and  body  model  was  varied  from  -8 
to  50  degrees.  The  lift,  drag,  Cj,,  and  Cp  were  then  plotted 
for  the  baseline  run.  Results  from  the  baseline  run  were  used 
to  determine  the  angles  of  attack  of  the  main  body  and  wing  to 
use  for  the  canard  runs.  Under  similar  tunnel  conditions  as 
the  baseline  run,  the  canard  was  varied  from  -25  to  +25 
degrees  deflection  for  various  wing/body  angles  of  attack  to 
determine  any  lift  enhancement  through  the  use  of  a  close- 
coupled  canard. 

B.  APPARATUS 

The  primary  equipment  used  was  the  NPS  low-speed  wind 
tunnel,  external  strain-gage  balance  and  signal  conditioning 
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assembly,  balance  calibration  rig,  canard/wing  model,  data 
acquisition  system,  and  data  reduction  software. 

1 .  Wind  Tunnel 

The  Naval  Postgraduate  School  low-speed,  single-return 
wind  tunnel  is  powered  by  a  100  hp  electric  motor  coupled  to 
a  three-blade  variable-pitch  fan  and  a  four-speed  truck 
transmission.  A  set  of  stator  blades  immediately  following 
the  fan  assist  in  straightening  the  flow.  A  combination  of 
turning  vanes  at  each  corner  and  two  fine-wire-mesh  screens  at 
the  entrance  to  the  settling  chairiber  help  reduce  the  air  flow 
turbulence  further.  The  settling  chamber  to  test  section 


contraction  ratio  is  about 

10:1.  [Ref. 

15] 

The  test  section 

of  the 

tunnel 

operates 

at 

approximately  atmospheric 

pressure 

due  to 

the  use 

of 

downstream  vents.  The  test  section  measures  45  by  32  inches. 
The  corner  lighting  and  reflection  plane  in  the  tunnel  test 
section  reduce  the  tunnel  height  from  32  to  28  inches,  which 
results  in  an  effective  cross-sectional  area  of  9.88  square 
feet.  A  remotely-controlled  turntable  mounted  flush  with  the 
reflection  plane  allowed  the  angle  of  attack  of  the  model  to 
be  varied.  The  temperature  of  the  tunnel  air  was  measured 
with  a  dial  thermometer  mounted  on  the  tunnel  wall  extending 
into  the  settling  chamber.  Figure  3  shows  the  NPS  low-speed 
wind  tunnel.  [Ref.  15] 
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Figure  3.  Naval  Postgraduate  School  Wind  Tunnel  [Ref.  15] 


The  test-section  dynamic  pressure,  q,  was  determined 
by  measuring  the  static  pressure  difference,  Ap,  between  the 
test  section  and  the  settling  chamber  using  a  water  manometer. 
The  settling  chamber  and  the  test  section  each  have  four 
static  pressure  taps  that  are  connected  to  the  manometer  via 
a  common  manifold.  The  pressure  difference  measured  by  the 
manometer,  in  centimeters  of  water,  was  converted  to  the  test- 
section  dynamic  pressure  and  test-section  reference  velocity 
from  a  previous  calibration  resulting  in  equations  (1)  and 
(2) .  [Ref.  16  and  Ref.  17] 

(1) 

g=  2.047  (-.026749  +  1.1149AP) 

(2) 

Where: 

p  air  density  (slugs/ft^3) 

Ap  manometer  reading  in  cm  of  H2O 

q  test-section  dynamic  pressure  (Ibf/ft^) 

2.047  a  constant  converting  cm  of  HjO  to  Ibf/ft^ 

1.1149  tunnel  calibration  factor 

-.026749  tunnel  calibration  intercept 

V  reference  velocity  (ft/sec) 

The  wind  tunnel  calibration  factor, 1 . 1149,  and  tunnel 
calibration  intercept,  -0.026749,  corrected  the  manometer 
reading,  Ap,  to  test -sect ion  dynamic  pressure,  q.  The 
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calibration  factor  was  found  by  plotting  the  actual  dynamic 
pressure  measured  by  a  pitot-static  tube  mounted  in  the  test 
section  versus  the  measured  pressure  difference.  The 
relationship  was  found  to  be  linear,  with  the  slope  of  the 
curve  being  the  tunnel  calibration  factor.  The  slope  did  not 
pass  through  the  origin,  which  resulted  in  there  being  a 
tunnel  calibration  intercept  with  the  y-axis.  [Ref.  17] 

2.  Canard/Ning  Model 

The  canard/wing  model  was  designed  as  a  half-model  for 
compatibility  with  the  existing  reflection-plane  balance 
previously  installed  in  the  wind  tunnel.  The  half-model  was 
of  a  generic  agile-fighter  fuselage  with  a  low-aspect-ratio 
close-coupled  canard  and  wing.  The  model  was  fabricated  from 
mahogany  by  Naval  Postgraduate  School  personnel.  The  canard 
and  wing  of  the  model  were  reinforced  with  an  aluminum  core. 
There  were  three  main  sections  to  the  model:  the  ogive  nose, 
canard,  and  wing.  The  ogive  nose  and  wing  section  were 
permanently  attached  to  one  another  by  an  aluminum  base  plate. 
A  large  removable  canard  section  allowed  the  canard  to  be 
instrumented  with  an  electric  motor  and  controller  so  the 
angle  of  attack  of  the  canard  could  be  varied  remotely  without 
shutting  down  the  tunnel.  The  model  mounted  flush  to  the  base 
of  the  tunnel  reflection  plane.  The  model  angle  of  attack  was 
varied  using  the  tunnel  turntable.  Figure  4  shows  a  sketch  of 


the  model.  See  Appendix  A  for  the  model  design  process  and 
the  resultant  geometric  parameters. 


3.  Balance  and  Turntable 

An  external  strain-gage  balance  and  turntable,  shown 
in  Figure  5,  was  originally  designed  and  built  by  NFS 
personnel  in  1974  to  facilitate  the  measurement  of  normal  and 
axial  forces  and  pitching  moment  in  the  NFS  low-speed  wind 
tunnel.  Each  external  strain-gage  bridge  had  four  active  legs 
for  automatic  temperature  compensation.  The  normal  and  axial 
moments  were  measured  by  two  orthogonal  strain-gage  bridges 
cemented  on  the  balance  colxiron  separated  by  a  vertical 
distance  of  26.5  inches.  With  the  wind  tunnel  in  operation 
the  force  on  the  model  created  a  different  moment  on  the  upper 
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and  lower  strain-gage  bridges.  Once  the  balance  was 
accurately  calibrated  the  voltage  output  from  the  lower  and 
upper  normal  or  axial  bridges  could  be  converted  to  moments 
and  subtracted  from  one  another,  then  divided  by  the  vertical 
separation  of  the  bridges,  to  find  the  normal  or  axial  force. 
Any  inherent  moments  on  the  model  were  removed  from  the 
measurements  by  this  process . 

Unfortunately,  prior  to  this  thesis,  the  balance  had 
never  been  calibrated  and  some  of  the  documentation,  such  as 
the  gage  factor  of  the  strain  gages,  was  unavailable. 
Drawings  of  the  balance  were  found  which  gave  the  separation 
distance  between  the  strain  gages  and  indicated  that  the 
balance  was  designed  to  be  capable  of  measuring  forces  up  to 
150  lbs.  A  balance  calibration  procedure  and  associated 
calibration  rig  were  then  developed  with  the  help  of  NPS 
personnel  and  Reference  18.  The  calibration  procedure  and 
associated  calibration  rig  are  described  in  Appendix  B. 

The  balance  and  turntable  were  one  rigid  unit.  When 
the  turntable  rotated  the  balance  column  upon  which  the  strain 
gages  were  cemented  rotated  with  it.  The  turntable  was 
controlled  by  an  electric  motor  with  hard-wired  remote.  The 
model  was  attached  to  the  top  of  the  turntable  platter.  The 
angle  of  attack  of  the  model  was  varied  using  the  tunnel 
turntable . 
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Figure  5.  Tunnel  Reflection  Plane  Balance  and  Turntable 


4.  Data  Acquisition  Hardwara 

Each  strain-gage  bridge  had  an  individual  signal 
conditioning  assembly  that  supplied  the  excitation  voltage. 
The  signal  conditioning  assemblies  allowed  their  associated 
strain-gage  bridges  to  be  zeroed  and  calibrated.  Each  balance 
channel  from  the  signal  conditioner  assembly  was  passed 
through  a  1000-gain  low-noise  amplifier.  The  signal  was  then 
routed  to  the  Hewlett-Packard  relay  multiplexer  that  sampled 
each  channel  sequentially  every  0.9  seconds.  The  sampling 
period  of  the  multiplexer  could  be  varied,  but  it  was  found 
that  a  sampling  period  less  than  0.9  seconds  caused  an 
unacceptable  level  of  noise  in  the  voltage  output.  A  sampling 
period  greater  than  0 . 9  seconds  did  not  measurably  reduce  the 
noise  in  the  voltage  output.  A  Hewlett-Packard  digital 
multimeter  then  converted  the  voltage  output  from  each  channel 
from  analog  to  digital.  An  IBM-AT  microcomputer  drove  the 
data  acquisition  software  and  Hewlett-Packard  hardware  and 
stored  collected  data.  A  sketch  of  the  data  acquisition 
system  is  shown  in  Figure  6.  [Ref.  15] 
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5.  Data  Acquiaition  Software 

The  data  acquisition  software  consisted  of  two  BASIC 
programs.  The  first  program,  JMK.BAS,  was  composed  of  a  shell 
program  to  control  the  Hewlett-Packard  instrumentation  and  a 
BASIC  program  to  manipulate  the  output  signal  from  the  data 
acquisition  hardware  into  a  usable  form.  The  voltage  outputs 
of  each  of  the  balance  channels  were  sampled  10  times  and  went 
through  a  standard  deviation  rejection  routine.  Any  reading 
that  fell  one  standard  deviation  outside  was  rejected.  The 
remaining  voltage  readings  for  each  channel  were  then 
averaged.  At  most  4  readings  were  rejected  of  the  10  sampled. 
In  general,  the  rejected  channel  readings  varied  from  the 
average  channel  readings  by  1  to  10%.  The  average  channel 
readings  were  used  to  calculate  the  normal  and  axial  force 
using  equations  found  in  Appendix  B.  The  normal  and  axial 
forces  at  the  given  angle  of  attack  were  used  to  calculate  the 
lift  and  drag  forces  using  equations  (3)  and  (4) . 

(3) 

DRAG  =  (Axial  Force)  sin  (AOA)  -  (Normal  Force)  cos  (AOA) 

(4) 

LIFT  =  (Axial  Force)  COS  (AOA)  +  (Normal  Force)  sin  (AOA) 

The  data  from  the  run  were  stored  on  the  C  drive  of 
the  IBM  microcomputer  and  on  a  floppy  disk  in  the  A  drive. 
Upon  completion  of  the  tunnel  run,  the  program  prompted  the 


user  to  input  the  test  conditions  for  the  tunnel  run  which 
were  added  to  the  files  on  the  C  and  A  drives  of  the 
microcomputer.  The  files  generated  by  JMK.BAS  were  later 
manipulated  by  the  program  CORCOEF.BAS  to  convert  lift  and 
drag  forces  into  Cl  and  Co  after  accounting  for  test  conditions 
and  making  necessary  corrections  for  tunnel  blockage,  balance 
calibration,  and  turntable  alignment.  The  data  acquisition  and 
manipulation  programs  are  listed  in  Appendices  C  and  D.  [Ref. 
15] 

C.  EXPERIMENTAL  CONDITIONS 

There  are  numerous  variables  which  could  affect  flow 
separation  and  vortex  formation  at  high  angles  of  attack.  The 
following  parameters  were  kept  as  constant  as  possible: 

(1)  Test  section  AP  *  17  cm  of  HjO 

(2)  Test  section  velocity  =  172  “174  ft/sec 

(3)  Test  section  Mach  number  =  0.17 

(4)  Reynolds  number  =  9.5  x  10®  based  on  wing  MAC 

A  strong  dependence  on  test-section  temperature  was  found 
for  one  of  the  balance  strain  gages.  The  application  of 
excessive  temperature  corrections  was  avoided  by  keeping  the 
tunnel  operating  temperature  below  74  degrees.  If  the  tunnel 
temperature  rose  above  this  value,  the  tunnel  was  shut  down 
and  allowed  to  cool.  See  the  Experimental  Corrections  section 
for  strain-gage  temperature  correction  procedures. 
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A  1/8-inch  gap  existed  between  the  base  plate  of  the  model 
and  the  reflection  plane.  The  gap  was  needed  to  prevent  the 
model  from  transferring  any  wind-induced  loads  to  the 
reflection  plane.  The  1/8-inch  gap  was  used  for  both  the 
wing/body  and  canard/wing/body  studies.  No  correction  was 
applied  for  the  gap;  however,  the  gap  should  not  adversely 
affect  the  experimental  results  since  the  experiment  was 
attempting  to  measure  the  comparative  lift  enhancement  between 
the  two  configurations. 

The  presence  of  the  wind  tunnel  walls  leads  to  stream-line 
curvature  and  downwash  errors  requiring  corrections  for  angle 
of  attack  and  induced  drag  of  the  model.  The  solid  boundaries 
form  streamlines  forced  upon  the  flow  which  result  in  the 
formation  of  image  vortices  which  should  be  accounted  for  in 
measurements  of  absolute  lift  and  drag.  In  the  present  study, 
the  errors  are  expected  to  be  small,  since  the  model  spans 
only  37.8%  of  the  test  section.  Since  the  study  involved  a 
comparison  of  the  lift  and  drag  for  the  wing/body  and 
canard/wing/body  configurations,  the  lack  of  application  of 
wall  corrections  should  not  adversely  impact  the  results. 

The  pitching  moment  of  the  model  was  not  measured.  The 
goal  of  this  thesis  was  to  measure  and  compare  the  enhanced 
lift  between  a  wing/body  and  canard/wing/body  configuration. 
As  this  study  was  of  a  basic  research  phenomenon,  no  attempt 
was  made  to  trim  the  aircraft  longitudinally.  A  strain-gage 
bridge  existed  that  when  calibrated  would  have  allowed  the 
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pitching  moment  of  the  model  to  be  measured.  If  in  the  future 
it  is  desired  to  measure  the  pitching  moment  of  the  model, 
then  the  moment  strain-gage  bridge  can  be  used. 

During  wind  tunnel  operation,  the  vibration  of  the  tunnel 
itself  could  not  be  controlled.  This  wind  tunnel  vibration 
was  transferred  to  the  tunnel  balance  via  the  model.  The 
electrical  outputs  of  the  strain-gage  bridges  were  averaged 
over  time  and  then  used  to  find  the  forces  on  the  model. 
Because  the  average  values  of  the  electrical  outputs  were 
used,  any  adverse  effect  due  to  vibration  of  the  tunnel  should 
have  been  ameliorated. 

The  ambient  turbulence  level  of  the  wind  tunnel  was  about 
0.2%.  The  reference  length  used  to  calculate  the  Reynolds 
number  was  the  mean  aerodynamic  chord  of  the  wing  taken  to  the 
centerline  of  the  fuselage  of  9.52  inches. 

D.  EXPERIMENTAL  PROCEDURE 

1.  Prerun  Calibration  and  Test 

The  external  strain-gage  balance  was  initially 
calibrated.  Appendix  B  illustrates  the  procedure  that  was 
used  to  find  the  coefficients  for  the  calibration  matrices. 
Once  the  balance  was  calibrated  it  was  checked  prior  to  each 
use . 

Figure  7  shows  a  sketch  of  the  calibration  rig 
assembly.  First  the  calibration  rig  was  attached  to  the 
balance  turntable  platter.  The  stand  and  cable  were  lined  up 
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vertically  and  horizontally  with  the  calibration  rig  head  by 
the  use  of  a  plumb  bob  and  level.  Proper  wiring  between  the 
signal  conditioning  assembly  and  the  HP  digital  multimeter  was 
checked.  Table  1  shows  the  channels  and  their  associated 
strain-gage  bridges.  The  IBM-AT  microcomputer  was  then 
energized  and  the  HP  program  PANELS.EXE  for  the  multimeter  was 
called.  Relay  Mux  1.  was  entered  in  the  PANELS.EXE  program 
and  channel  (8)  was  enabled.  The  Pacific  Amplifier  gain 
switch  was  turned  to  zero  and  the  gain  output  was  adjusted  via 
a  set  screw  to  +/-  50  jivolts.  The  gain  was  increased  to  1000 
and  the  gain  input  adjusted  to  +/-  500  ^ivolts.  Channels  (2) 
through  (5)  were  read  and  recorded.  In  general,  the  voltage 
readings  from  channels  (2)  to  (5)  were  not  zero.  The  signal 
conditioner  assemblies  could  have  been  adjusted  so  that  the 
outputs  from  the  channels  were  close  to  zero.  Adjusting  the 
voltage  readings  of  the  channels  risked  unsettling  the  system 
or  changing  the  calibration  of  the  balance.  It  was  decided  to 
note  the  voltage  readings  from  channels  (2)  through  (5)  as 
zero  offsets. 

Initially  no  weight  was  attached  to  the  calibration 
rig  assembly.  The  IBM  BASIC  editor  was  entered  and  the 
program  JMK.BAS  called.  The  program  prompted  the  user  for  the 
angle  of  attack  of  the  model,  the  deflection  angle  of  the 
canard,  the  temperature,  and  the  zero  offset  readings  of 
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channels  (2)  to  (5)  found  using  PANELS.EXE.  The  requested 
information  was  entered  and  the  program  started.  The 
displayed  axial  and  normal  forces  found  should  be  less  than 
.05  Ibf.  If  the  resultant  normal  and  axial  forces  were 
greater  than  .05  Ibf,  then  the  offset  voltages  from  channels 
(2)  through  (5)  were  checked  and  reentered.  It  was  not 
necessary  to  leave  the  program  JMK.BAS  and  reenter  PANELS.EXE 
to  update  the  zero  offset  voltages  because  the  program 
displayed  the  average  voltage  readings  for  each  channel  on  the 
monitor.  These  displayed  readings,  that  with  no  force  on  the 
calibration  rig  were  the  zero  offsets,  were  then  entered  into 
the  computer  at  the  start  of  the  next  iteration. 

Once  the  balance  was  zeroed  the  turntable  was  rotated 
to  0,  90,  and  55  degrees.  Suspending  weights  from  the 
calibration  rig  with  the  turntable  at  0  degrees  imparted  a 
pure  normal  moment  to  the  balance.  Rotating  the  turntable  to 
90  degrees  and  suspending  weights  induced  a  pure  axial  moment 
on  the  balance.  The  55-degree  position  was  used  to  impart  a 
simultaneous  axial  and  normal  moment  on  the  balance.  The  55- 
degree  position  simulated  the  approximate  resultant  force 
direction  on  the  balance,  which  was  at  about  325  degrees, 
while  the  tunnel  was  in  operation  with  the  model  at  high 
angles  of  attack.  With  the  turntable  rotated  to  one  of  the 
angles  of  attack  the  requested  information  was  entered  into 
the  computer  including  the  previously  determined  zero  offset 
for  each  channel.  Successively  larger  weights  were  then  hung 
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from  the  balance  calibration  rig  and  the  normal  and  axial 
forces  calculated  by  the  program  JMK.BAS.  After  each  weight 
was  hung  from  the  balance  the  zero  offsets  of  the  channels 
were  checked  to  ensure  that  they  had  not  drifted.  In  general, 
channel  (3)  drifted  1  or  2  millivolts  a  minute.  If  the  zero 
offset  for  channel  (3)  was  not  updated  over  the  course  of  one 
hour,  error  in  the  balance  would  change  by  2%.  The  computed 
values  for  axial  and  normal  forces  were  then  compared  to  the 
actual  normal  and  axial  forces  induced  by  the  weight  on  the 
balance  rig  and  the  computed  values  found  during  the  original 
calibration  procedure.  (See  Appendix  B.)  If  the  computed 
values  for  the  prerun  calibration  differed  from  the  original 
computed  calibration  values  by  less  than  .1  Ibf,  alignment  of 
the  turntable  was  verified  and  a  known  weight  was  again 
suspended  from  the  rig.  If  the  difference  was  still  greater 
than  .1  Ibf,  the  balance  was  recalibrated. 

The  calibration  rig  was  sized  for  a  maximum 
anticipated  force  on  the  balance  of  about  25  Ibf  during  tunnel 
operation.  It  was  found  that  the  balance  was  actually 
subjected  to  axial  forces  of  up  to  70  Ibf.  Due  to  the 
limitations  of  the  calibration  rig  the  largest  weight  that 
could  be  hung  from  the  rig  was  21  Ibf.  The  response  of  the 
strain  gages  was  linear,  so  the  use  of  lesser  weights  should 
still  verify  the  calibration. 
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2 .  Data  Collection 


After  a  satisfactory  prerun  calibration  the  model  was 
mounted  on  the  turntable  platter.  If  the  canard  was  used, 
the  canard  remote-control  wiring  was  strung  through  the 
turntable  platter  and  balance  column.  The  program  PANELS.EXE 
was  then  called  and  the  Pacific  Amplifier  was  again  zeroed  in 
the  same  manner  as  was  done  during  the  prerun  calibration. 
The  outputs  of  channels  (2)  to  (5)  were  recorded  for  zero 
offsets . 

The  program  JMK.BAS  was  called  via  the  BASIC  editor 
and  the  requested  information  was  entered.  The  zero  offsets 
for  channels  (2)  to  (5)  were  then  updated  by  letting  the 
program  run  before  the  tunnel  was  started.  The  only  offset 
that  usually  changed  by  a  significant  amount  was  for  channel 
(3) .  Tunnel  runs  were  made  changing  the  turntable  angle  of 
attack  and/or  canard  deflection  angle  as  required.  The  remote 
control  of  the  canard  allowed  the  canard  deflection  angle  to 
be  varied  without  actually  shutting  down  the  tunnel  to 
reposition  the  canard. 

3 .  Preliminary  Runs 

Test  runs  were  made  at  a  number  of  different  tunnel 
speeds  with  the  wing/body  model  to  verify  the  operation  of  the 
balance.  A  number  of  short  runs  were  made  during  testing  of 
the  canard  servo  motor  and  eventual  canard  electric  motor  with 
lead  screw.  As  noted  in  Appendix  A,  the  servo  motor  was 
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unable  to  position  and  hold  the  canard  in  the  flow  at  all  but 
the  lowest  tunnel  speeds,  which  necessitated  the  use  of  an 
electric  motor  with  lead  screw  to  drive  the  canard. 

It  was  found  that  there  was  a  linear  dependence  of 
channel  (3)  upon  temperature  for  temperatures  greater  than 
about  74  degrees-  In  general,  no  data  were  used  that  were 
taken  at  a  tunnel  settling  chamber  temperature  greater  than  74 
degrees.  Once  the  tunnel  reached  this  temperature  the  tunnel 
was  shut  down  to  allow  it  to  return  to  ambient  air 
temperature.  At  angles  of  attack  of  40  and  48  degrees,  the 
tunnel  heated  up  very  rapidly  from  the  typical  starting 
temperature  of  67  degrees.  The  time  until  the  tunnel  reached 
74  degrees  was  as  little  as  15  minutes.  The  relative  slowness 
of  the  data  acquisition  system  meant  as  few  as  six  data  points 
could  be  collected  under  such  conditions.  Accordingly,  a 
compensation  routine  was  incorporated  in  the  program  JMK.BAS 
to  correct  for  the  temperature  dependence  and  allow  for  longer 
tunnel  runs  should  the  need  arise  in  the  future.  Table  2  and 
Figure  8  show  the  linear  variation  of  the  channel  (3)  voltage 
output  with  temperature  at  a  constant  model  angle  of  attack  of 
40  degrees,  a  constant  canard  deflection  angle  of  -17  degrees, 
and  a  AP  of  17  cm  H2O. 
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E .  EXPERIMENTAL  CORRECTIONS 

The  lift  and  drag  coefficients  were  corrected  for  model 
blockage  by  equations  (5)  and  (6)  from  References  15  and  18. 

(5) 

q  =  q„(l  +  2e) 

U  =  U„(l  +  e)  (6) 

where : 

q  dynamic  pressure  (Ibf/ft^) 

q„  measured  reference  q  (Ibf/ft*) 

U  horizontal  velocity  (ft/sec) 

Um  measured  horizontal  velocity  (ft/sec) 

e  blockage  factor 

The  blockage  correction  factor  depended  upon  the  angle  of 
attack  of  the  model  and,  if  the  canard  was  in  .^lace,  the 
canard  deflection  angle.  For  the  wing/body  case^  the  blockage 
correction  factor  was  determined  from  e*.  nation  (7)  where 
0.02098  was  the  ratio  of  the  cross-sections  of  the  wing/body 
and  the  tunnel  test  section.  If  the  canard  was  in  place, 
equation  (8)  was  used  where  0.00383  was  the  ratio  of  the 
canard  cross-sectional  area  to  that  of  the  test  section. 

Efcody  *  •0298sin(AOA)  (7) 
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(8) 


Ewycanard  =  ^body  +  •  00383sin  (AOA  +  Canard  AOA) 

During  tunnel  postrun  calibration  trials,  a  calibration 
error  was  found  in  the  balance.  The  balance  calibration  error 
was  quantified  using  the  calibration  rig  with  the  turntable  at 
an  angle  of  attack  of  55  degrees.  With  the  turntable  at  55 
degrees  a  resultant  force  direction  of  325  degrees  could  be 
simulated,  which  was  close  to  the  resultant  force  direction  on 
the  model  at  high  angles  of  attack.  Figure  9  is  a  plot  of  the 
experimental  normal  force  versus  actual  normal  force  at  an  AOA 
of  55  degrees.  Figure  10  is  a  plot  of  the  experimental  axial 
force  versus  the  actual  axial  force  at  an  AOA  of  55  degrees. 
The  plots  indicate  an  error  of  +4.5%  in  the  reading  of  the 
actual  normal  force  and  an  error  of  -3.7%  in  the  reading  of 
the  actual  axial  force.  Accordingly,  the  program  CORCOEF.BAS 
was  updated  to  multiply  the  normal  force  output  by  . 955  and 
the  axial  force  output  by  1.037  to  compensate  for  balance 
calibration  error  these  values  were  then  used  to  find  the 
corrected  lift  and  drag.  Table  3  includes  the  data  used  in 
Figures  9  and  10. 

The  turntable  was  found  to  be  misaligned  by  +2.2  degrees 
at  the  end  of  all  the  data  runs.  This  necessitated  a  further 
correction  of  the  calculated  lift  and  drag  forces.  The  model 
was  mounted  on  the  platter  with  the  turntable  at  90  degrees 
AOA  with  the  result  that  a  -2.2  degree  correction  was  applied 
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to  the  AOA  of  the  wing/body.  All  the  corrections  were 
incorporated  in  the  BASIC  program  CORCOEF.BAS  found  in 
Appendix  D. 
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Figure  9.  Experimental  vs.  Actual  Normal  Force 
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Table  3.  Experimental  and  Actual  Axial  and  Normal  Forces 
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10.2216 


13.2216 


14.2262 


15.2251 


16.2297 


17.2241 


18.2287 


21.2388 


6.11 


8.02 


8.48 


9.16 


9.74 


10.30 


10.99 


11.40 


5.86 


7.58 


8.16 


8.73 


9.31 


9.88 


10.46 
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AXIAL  EXPER 
FORCE  (LBS) 


8.03 


10.50 


11.16 


12.04 


12.83 


13.54 


14.44 


15.09 


16.81 


8.37 


10.83 


11.65 


12.47 


13.29 


14.11 


14.93 


15.75 
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III.  RESULTS  AND  DISCUSSION 


The  following  sections  discuss  the  plots  and  data  from  six 
wind-tunnel  runs.  The  first  was  a  baseline  run  without  the 
canard,  where  the  angle  of  attack  (AOA)  of  the  wing/body  was 
varied  from  -8  to  50  degrees  in  2-degree  increments  or  less. 
The  following  runs  were  conducted  with  the  canard  in  place  at 
wing/body  angles  of  attack  of  10,  22,  34,  40,  and  48  degrees 
where  the  canard  deflection  angle  (canard  AOA)  was  varied  from 
-25  to  25  degrees  in  increments  of  5  degrees  or  less.  It 
should  be  noted  that  the  term  "canard  deflection  angle"  refers 
to  the  incidence  angle  of  the  canard  relative  to  the  fuselage 
centerline.  Ci  and  Co  refer  to  the  coefficients  of  lift  and 
drag.  Data  for  Cj,  versus  AOA,  Cj,  versus  Co,  and  Cl/Cd  versus 
AOA  were  plotted  for  the  baseline  run.  Data  for  Cj,  versus 
canard  deflection  angle  and  Ci,  versus  were  plotted  for  the 
five  wing/body  angles  of  attack  with  the  canard  in  place. 

The  reference  area  for  the  Cj,  and  the  Cp  of  the  wing/body 
configuration  was  the  area  of  the  wing  taken  to  the  model 
centerline  of  0.685  ft^.  The  reference  area  for  the  Cj,  and  the 
Cd  of  the  canard/wing/body  configuration  was  the  area  of  the 
wing  to  the  centerline  plus  the  exposed  area  of  the  canard, 
which  total  area  equalled  0.815  ft*.  Appendix  E  includes  the 
data  sets  used  for  all  plots  in  the  Results/Discussion 
chapter . 
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A.  BASELINE  WING/BODY  CONFIGDBATION 

Figure  11  shows  the  variation  of  Cl  with  AOA.  The  AOA  of 
the  wing/body  was  varied  from  -8  to  50  degrees, and  Cl  varied 
from  .036  to  1.586.  The  maximum  lift  occurred  at  an  AOA  of  40 
degrees . 

Wings  with  aspect  ratios  less  than  5  typically  stall  at 
higher  angles  of  attack  than  higher-aspect -ratio  wings.  The 
wing  of  the  model  had  an  aspect  ratio  of  3.  The  overall  shape 
of  the  Cl  versus  AOA  curve  is  similar  to  curves  for  other  low- 
aspect-ratio  wing  configurations.  [Ref.  19] 

The  Cl  versus  AOA  curve  was  linear  from  -8  to  4  degrees 
with  a  slope  of  0.0577/degree.  At  angles  of  attack  greater 
than  4  degrees  the  slope  became  slightly  non-linear;  from  4 
to  12  degrees,  the  Cl  was  greater  than  the  projected  linear 
slope.  This  non-linear  increase  in  lift  was  probably  due  to 
the  formation  of  a  weak  wing-leading-edge  vortex.  From  12  to 
18  degrees,  the  Cl  was  less  than  the  projected  linear  slope. 
At  about  18  degrees  the  first  maximum  Cl  occurred  and  the 
first  stall  began. 

The  angle-of-attack  increment  was  decreased  to  one 
degree  from  17  to  29  degrees  so  that  the  nature  of  this  first 
stall  could  be  better  observed.  It  was  established  that  from 
21  to  25  degrees  the  Cl  did  not  appreciably  vary  with  AOA. 
The  first  stall  was  probably  caused  by  the  onset  of  major  flow 
separation  on  the  outboard  wing  section.  At  angles  of  attack 
greater  than  25  degrees,  the  Cl  again  began  to  increase  with 
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AOA.  It  is  speculated  that  a  strong  leading-edge  vortex  has 
formed  at  this  high  angle  of  attack,  tending  to  reattach  the 
flow. 

From  28  to  32  degrees,  the  Cj,  versus  AOA  curve  was  again 
linear  with  a  slope  of  0.0473/degree.  This  slope  is  not  very 
different  from  the  earlier  linear  slope  from  -8  to  4  degrees. 
At  angles  of  attack  greater  than  about  32  degrees  the 
versus  AOA  curve  again  began  to  decrease  as  the  second  stall 
was  initiated.  The  maximum  Ci,  was  1.586  at  an  AOA  of  40 
degrees.  The  0^  at  40  degrees  was  0.53  greater  than  the 
during  the  first  stall.  The  second  stall  was  probably  due  to 
the  breakdown  of  the  strong  wing-leading-edge  vortex. 

Figure  12  shows  the  variation  of  Cl  with  C^.  The  value 
of  Co  at  the  maximum  Cl  was  1.262.  The  value  for  C^  minimum 
was  0.034  at  a  Cl  of  0.036.  The  plot  shows  typical  drag-polar 
behavior  from  a  Cl  of  -0.420  to  1.031.  From  a  Cl  of  1.031  to 
1.059,  the  first  stall  occurred  and  Cp  increased  from  .390  to 
.512.  For  a  Cl  greater  than  1.059,  lift  again  increased  until 
the  maximum  Cl^  1.586,  was  reached.  At  a  Cp  greater  than  1.262 
the  second  stall  took  place  where  the  lift  fell  off  rapidly 
and  the  drag  increased  greatly. 

Figure  13  show  the  variation  of  Cl/C^  with  AOA.  The 
maximum  Cl/Cd  was  7.15  at  an  AOA  of  6  degrees.  This  condition 
is  the  region  in  which  the  lift  was  enhanced  by  the  weak 
leading-edge  vortex.  The  minimum  post-stall  Cl/Cd  was  0.92  at 
50  degrees.  Cl/Cq  increased  rapidly  from  -8  to  6  degrees. 
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changing  12.75  units  in  14  degrees.  From  6  degrees  to  50 
degrees  C^/Cb  decreased  in  a  smooth  non-linear  drop-off  lo  0.92 
at  an  AOA  of  50  degrees.  Values  in  the  high  angle-of-attack 
regime  from  20  to  40  degrees,  where  an  agile-aircraft  might  be 
expected  to  have  an  excursion,  decrease  from  about  2.8  to  1.3, 
indicating  the  significant  thrust  levels  required  for 
supermaneuverability . 
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Figure  12.  Lift  vs  Drag  Coefficient;  Baseline 


Figure  13.  Ci/Co  vs.  Angle  of  Attack;  Baseline 
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B.  CANARD/WING/BODY  CONFIGUBATION;  AOA  10  DEGREES 

Figure  14  shows  the  variation  of  Ci,  with  canard  deflection 
angle  at  a  wing/body  AOA  of  10  degrees.  The  maximum  Cl  was 
0.698  at  17  degrees.  From  -25  to  0  degrees  the  increase  in  Cl 
with  canard  deflection  angle  was  fairly  linear.  Beyond  0 
degrees  the  curve  was  non-linear  with  an  increased  data 
scatter.  Lift  began  to  decrease  at  deflection  angles  greater 
than  20  degrees. 

When  compared  to  the  baseline  wing/body  configuration, 
lift  was  slightly  enhanced  for  a  canard  deflection  angles 
between  7  and  about  20  degrees.  Positioning  the  canard 
outside  of  this  range  actually  caused  a  lift  degradation. 
Lift  was  maximized  at  a  canard  angle  of  17  degrees,  where  the 
maximum  Cl  was  3.41%  greater  than  the  baseline  value  of  0.675. 
At  a  canard  deflection  angle  of  25  degrees  Cl  was  1.48%  less 
than  the  baseline  value.  It  can  be  seen  that  the  lift 
enhancement  with  the  addition  of  a  canard  is  slight  at  the 
relatively  low  wing/body  AOA  of  10  degrees. 

Figure  15  shows  the  variation  of  Cl  with  C^  for  this  case. 
The  maximum  Cl  was  0.698  at  a  C^  of  0.188.  Co  minimum  was 
0.094  at  a  Cl  of  0.631.  At  a  Cq  of  0.120  and  0.143,  Cl  dropped 
slightly  compared  to  the  data  trend.  These  two  data  points 
were  repeatable  and  were  not  due  to  scatter.  At  a  Co  of  about 
0.2  lift  began  to  drop  off  markedly  and  drag  increased 
greatly. 
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C.  CANARD /HING/BODY  CONFIGURATION;  AOA  22  DEGREES 

Figure  16  shows  the  variation  of  Cj.  with  canard  deflection 
angle  for  a  wing/body  AOA  of  22  degrees.  The  maximum  was 
1.422  at  a  canard  angle  of  7  degrees.  There  was  little  data 
scatter  and  the  curve  was  fairly  smooth  for  all  canard 
deflection  angles  tested.  The  lift  varied  in  a  strong  non¬ 
linear  fashion  from  a  Cj.  of  0.97  to  1.36.  A  typical  stall 
behavior  was  indicated,  the  lift  coefficient  leveling  off  to 
1.30-1.35  after  stall. 

When  compared  to  the  baseline  configuration,  lift  was 
greatly  enhanced  for  a  canard  deflection  angle  between  -15  and 
25  degrees.  Positioning  the  canard  from  -25  to  about  -15 
degrees  degraded  the  lift  below  the  baseline  value.  As  the 
canard  angle  was  increased  from  -15  to  7  degrees,  lift  was 
enhanced  significantly.  Lift  was  maximized  at  a  canard  angle 
of  7  degrees,  where  the  Cl  was  34%  greater  than  the  baseline 
value  of  1.061.  From  7  to  25  degrees  the  values  of  Cl 
decreased  with  an  increase  in  canard  angle;  but,  the 
comparative  lift  enhancement  was  still  significant.  At  a 
canard  deflection  angle  of  25  degrees.  Cl  was  still  24.5% 
greater  than  the  baseline  value.  It  should  be  noted  that  this 
enhancement  is  taking  place  at  the  wing/body  angle  of  attack 
where  the  first  stall  occurred.  Evidently  the  canard/wing 
interaction  is  preventing  the  wing  separation  from  taking 
place,  the  canard  vortex  providing  the  energizing  mechanism. 
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Figure  17  shows  the  variation  of  Cj,  with  Cd.  The  maximum 
Cj,  was  1.422  at  a  Co  of  0.547.  Cq  minimum  was  0.368  at  a  of 
1.007.  The  overall  plot  was  relatively  smooth  with  little 
data  scatter.  The  slope  was  approximately  linear  from  a  Co  of 
0.413  to  0.502.  At  a  Cj,  greater  than  1.422  stall  occurred  and 
lift  dropped  of  markedly  while  drag  increased. 


LIFT  COEFFICIENT 


Figure  16.  Lift  Coefficient  vs.  Canard  Deflection; 
AOA  22 
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D.  CANARD /WIN6/B0DY  CONFIGURATION;  AOA  34  DEGREES 

Figure  18  shows  the  variation  of  Cj,  with  canard  deflection 
angle  for  a  wing/body  AOA  of  34  degrees.  This  angle  of  attack 
is  in  the  region  midway  along  the  second  rise  of  the  baseline 
configuration  lift  curve.  The  maximum  Ci.  was  1.642  at  a 
canard  angle  of  -7  degrees.  The  Cj,  increased  with  canard 
angle  from  -25  to  -7  degrees.  At  -7  degrees  the  primary  Cj, 
peak  occurred;  Cl  then  decreased  with  increasing  canard  angle 
until  a  canard  deflection  angle  of  0  degrees.  The  then 
increased  with  increasing  canard  deflection  angle  until  an 
angle  of  5  degrees  was  reached.  At  5  degrees  there  was  a 
secondary  Cj,  peak  of  1.625.  At  canard  deflection  angles 
greater  than  5  degrees  the  canard/wing/body  began  to  stall. 
Within  a  band  of  data  scatter,  the  data  were  repeatable,  in 
particular  the  twin-peaked  behavior. 

From  -25  to  -10  degrees  lift  was  degraded  using  the 
canard.  As  the  canard  deflection  angle  was  increased  from  -12 
to  -7  degrees,  the  lift  was  enhanced  over  the  baseline  case. 
Lift  was  maximized  at  a  canard  angle  of  -7  degrees,  where  Ci, 
was  9.39%  greater  than  the  baseline  value  of  1.501.  From  -7 
to  0  degrees  Cl  decreased  with  canard  deflection;  but  the  lift 
was  still  greater  than  the  baseline  configuration.  Lift 
increased  with  increasing  canard  angle  from  0  to  5  degrees. 
At  a  canard  angle  of  25  degrees  Cl  was  2.53%  greater  than  the 
baseline  value. 
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Figure  19  shows  the  variation  of  Cj,  with  Cd.  The  maximum 
Cj,  was  1.642  at  a  Co  of  0.982.  The  plot  was  relatively  smooth 


and  linear  for  coefficients  of  drag  less  than  0.982.  At  a  Co 
of  0.982  the  data  became  multi-valued,  but  all  data  points 
were  shown  to  be  repeatable.  At  a  greater  than  1.058 
stall  occurred  and  lift  decreased  markedly  while  drag 
increased. 
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E.  CAMABD/WIN6/B0DY  CONFIGOBATION;  AOA  40  DEGREES 

Figure  20  shows  the  variation  of  Cj,  with  canard  deflection 
angle  for  a  wing/body  AOA  of  40  degrees.  This  wing/body  angle 
corresponds  to  the  condition  of  the  maximum  Cj,  for  the 
baseline  configuration.  The  maximum  Cj,  was  1.700  at  a  canard 
angle  of  -15  degrees.  The  Cl  increased  with  canard  deflection 
angle  and  the  slope  was  fairly  steep  from  -25  to  -15  degrees. 
At  canard  deflection  angles  greater  than  -15  degrees  the 
canard/wing/body  began  to  stall. 

Lift  was  enhanced  compared  to  the  baseline  configuration 
for  canard  deflections  between  -25  and  15  degrees.  Lift  was 
degraded  for  canard  angles  greater  than  15  degrees.  A  canard 
deflection  between  -20  and  5  degrees  resulted  in  at  least  a 
3.97%  increase  in  lift  over  the  baseline  value  of  1.586.  Lift 
was  maximized  at  a  canard  deflection  of  -15  degrees  ,  where  the 
Cl  was  7.19%  greater  than  the  baseline  value.  From  -15  to  0 
degrees  Cl  decreased  slightly  with  canard  deflection  to  1.680; 
but  overall  the  lift  was  still  greater  than  the  baseline.  Cl 
was  15.45%  less  than  the  baseline  value  at  a  canard  deflection 
of  25  degrees. 

Figure  21  shows  the  variation  of  Cl  with  Co-  The  maximum 
Cl  was  1.700  at  a  Cq  of  1.214.  Q,  minimum  was  1.099  at  a  Cl  of 
1.586,  though  this  minimum  is  due  to  the  limited  data  points 
taken  rather  than  to  an  infinite  slope.  The  plot  was 
relatively  linear  for  coefficients  of  drag  less  than  1.214. 
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At  a  Cd  of  1.240  the  data  became  slightly  scattered.  At  a  Cp 
of  1.240  the  canard/wing/body  began  to  stall  and  at  a  Cd  of 
1.325,  Cl  lift  dropped  precipitously  while  there  was  little 
change  in  C^. 
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Figure  21.  Lift  vs.  Drag  Coefficient;  AOA  40 


F.  CANARD/WING/BODY  CONFIGURATION;  AOA  48  DEGREES 

Figure  22  shows  the  variation  of  Ci,  with  canard  deflection 
angle  for  a  wing/body  AOA  of  48  degrees.  This  wing/body  angle 
is  deep  in  the  second  post-stall  region  for  the  baseline 
configuration.  The  maximum  0^  was  1.649  at  -17  degrees.  The 
shape  of  the  plot  was  similar  to  the  plot  shown  in  Figure  20 
for  the  wing/body  at  an  AOA  of  40  degrees.  Cl  increased  with 
canard  deflection  for  angles  less  than  -17  degrees.  At  canard 
deflection  angles  greater  than  -17  degrees,  stall  occurred  and 
the  Cl  decreased  rapidly  with  increased  canard  deflection. 

Lift  was  enhanced  compared  to  the  baseline  configuration 
for  canard  deflections  between  -25  and  20  degrees.  Lift  was 
degraded  for  canard  deflection  angles  greater  than  20  degrees. 
A  canard  deflection  between  -25  and  0  degrees  resulted  in  at 
least  a  13.92%  increase  in  lift  over  the  baseline  value  of 
1.394.  Lift  was  maximized  at  a  canard  deflection  of  -17 
degrees,  where  the  Cl  was  18.29%  greater  than  the  baseline 
value.  Cl  was  1.15%  less  than  the  baseline  value  at  a  canard 
deflection  of  25  degrees.  In  this  post-stall  regime,  the 
canard/wing  interaction  remains  effective  over  a  much  wider 
angle-of-attack  range  than  for  the  previous  cases. 

Figure  23  shows  the  variation  of  Cl  with  C^.  The  maximum 
Cl  was  1.64  9  at  a  Co  of  1.534.  For  coefficients  of  drag  less 
than  1.541,  Cl  varied  by  less  than  +/-  0.9%  from  1.635.  At  a 
Cd  of  1.623  the  canard/wing/body  began  to  stall  and  Cl 
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decreased  greatly  to  a  value  of  1.378  while  Cj,  actually 
decreased  by  0.05  to  1.573.  It  is  interesting  to  observe  that 
at  canard  deflection  angles  greater  than  0  degrees,  the  drag 
decreases  with  decreasing  lift,  an  anomaly  in  a  typical  post¬ 
stall  regime. 
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Figure  22.  Lift  Coefficient  vs.  Canard  Deflection; 
AOA  48 
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6.  CC»«PARISON  BETWEEN  WING/BODY  AND  CANABD/WING/BODY 

Figure  24  shows  the  comparative  lift  of  the  wing/body  and 
the  canard/wing/body  configuration.  As  predicted  by  Lacey, 
there  was  little  lift  enhancement  using  a  close-coupled  canard 
at  wing/body  angles  of  attack  less  than  18  degrees.  The 
maximum  Cj,  at  10  degrees  AOA  was  only  3.4%  greater  than  the 
baseline  value.  At  22  degrees  AOA,  where  the  first  stall 
occurred,  the  lift  was  greater  than  the  baseline  value  by  34%. 
This  dramatic  increase  in  lift  is  thought  to  be  due  to  the 
downwash  of  the  close-coupled  canard  impinging  upon  the 
flowfield  of  the  wing,  thereby  delaying  the  onset  of  flow 
separation.  At  34  degrees  AOA,  the  lift  was  9.4%  greater 
than  the  baseline  value.  At  40  degrees  AOA,  where  the  lift 
was  previously  maximized  for  the  wing/body  configuration,  the 
lift  was  7.2%  greater  than  the  baseline  value.  In  general, 
the  lift  enhancement  was  not  as  great  where  major  separation 
did  not  previously  exist  for  angles  of  attack  tested  in  the 
baseline  configuration.  At  48  degrees  AOA,  where  the  second 
stall  occurred,  the  lift  was  greater  than  the  baseline  value 
by  18.3%  This  great  increase  in  lift  is  thought  to  be  due  to 
the  vortex  of  the  close-coupled  canard  constructively 
interfering  with  the  vortex  of  the  wing,  thus  delaying  the 
onset  of  flow  separation. 
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Figure  24.  Lift  of  Wing/Body  and  Maximum  Lift  of 
Canard/Wing/Body  vs.  Angle  of  Attack 
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Figure  25  compares  the  lift  of  the  wing/body  to  the 
minimum  lift  of  the  canard/wing/body  and  to  the  lift  of  the 
canard/wing/body  with  the  canard  at  a  canard  fixed  deflection 
angle  of  -5  degrees.  To  illustrate  that  the  canard  must  be 
positioned  properly  for  each  wing/body  AOA,  the  minimum  lift 
obtained  when  the  canard  deflection  angle  was  varied  was  also 
plotted.  Not  too  surprisingly,  if  the  canard  was  positioned 
to  minimize  the  overall  lift,  the  lift  of  the  canard/wing/body 
was  less  than  that  of  the  wing/body  configuration.  But  at  an 
AOA  of  48  degrees,  with  the  canard  positioned  to  minimize 
overall  lift,  the  lift  was  only  1.15%  less  than  the  wing/body 
configuration.  In  fact,  at  an  AOA  of  48  degrees,  if  the 
canard  was  at  any  canard  deflection  angle  other  than  25 
degrees  the  lift  was  improved  over  the  baseline  value. 

Experimental  results  indicated  that  a  canard  deflection 
angle  could  be  chosen  that  enhanced  lift  at  a  number  of 
different  wing/body  angles  of  attack.  In  fact,  if  the  canard 
was  positioned  at  the  approximate  average  canard  deflection 
angle  that  maximized  lift  for  all  runs,  a  -5  degree  incidence 
angle  with  respect  to  the  centerline  of  the  fuselage,  the  lift 
was  enhanced  for  all  tested  wing/body  angles  of  attack  except 
at  10  degrees.  At  an  AOA  of  10  degrees,  the  lift  was  6.5% 
less  than  the  baseline  value.  The  lift  was  14.6%  greater  than 
the  baseline  value  at  an  AOA  of  22  degrees.  At  AOA's  of  34 
and  40  degrees,  the  lift  was  7.7%  and  6.6%  greater  than  the 
corresponding  baseline  values.  The  lift  was  15.5%  greater 
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than  the  baseline  value  at  an  AOA  of  48  degrees.  Positioning 
the  canard  at  -5  degrees  did  not  maximize  the  lift  for  any 
wing/body  AOA;  but,  if  it  is  desired  to  minimize  the  weight 
and  complexity  of  a  moveable  canard  system  while  increasing 
the  lift  over  a  large  AOA  range,  a  fixed  -5  degree  canard 
incidence  angle  is  a  good  potential  choice  for  the 
configuration  studied. 
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Figure  26  compares  the  baseline  drag  polar  for  the 
wing/body  configuration  with  the  drag  polar  for  the 
canard/wing/body  configuration  using  a  canard  deflection  that 
maximized  the  lift  for  each  wing/body  AOA.  For  all  tested 
angles  of  attack  except  10  degrees,  the  lift  was  greater  for 
the  canard/wing/body  than  for  the  wing/body  configuration  at 
the  same  drag  coefficient. 

Figure  27  compares  the  lift-to-drag  ratio  for  the 
wing/body  configuration  with  that  for  the  canard/wing/body 
configuration.  For  the  canard/wing/body,  the  Cl/Cd  at  each 
tested  AOA  was  greater  than  the  baseline  value  for  all  tested 
wing/body  angles  of  attack  except  at  10  degrees.  At  an  AOA  of 
10  degrees,  the  Cl/Cd  ratio  was  34.5%  less  than  the  baseline 
value.  The  Cl/Co  ratios  were  7.4%,  8.9%,  11.4%,  and  10.4% 
greater  than  the  baseline  values  at  angles  of  attack  of  22, 
34,  40,  and  48  degrees.  Essentially,  using  a  properly  located 
close-coupled  canard  exacted  a  penalty  in  only  at  low 
angles  of  attack,  below  a  lift  coefficient  of  0.35.  At  higher 
angles  of  attack,  the  increase  in  over  the  baseline  value 
was  as  great  as  71.4%. 


63 


Figure  26.  Lift  vs.  Drag  Coefficient  for  Wing/Body 
and  Canard/Wing/Body 
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Figure  27.  Ci/Co  vs.  AOA  for  Wing/Body  and 
Canard/Wing/Body 
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Figure  28  shows  the  variation  of  the  canard  absolute  AOA 
for  maximum  lift  of  the  canard/wing/body  configuration  with 
the  wing/body  AOA.  The  canard  absolute  AOA  was  measured  with 
respect  to  the  freestream  velocity,  as  opposed  to  the  fuselage 
reference  line.  The  average  canard  absolute  AOA  for  wing/body 
angles  of  attack  of  10,  22,  34,  and  40  degrees  was  26.8 


Figure  28.  Absolute  Canard  AOA  (deg)  vs.  Body  AOA  (deg) 

degrees.  The  canard  absolute  AOA  varied  only  by  +/-  2  degrees 
across  the  wing/body  angle-of-attack  range.  At  a  wing/body 
AOA  of  48  degrees,  the  canard  absolute  AOA  was  30.8  degrees, 
4  degrees  greater  than  the  previous  average.  The  Cj,  versus 
canard  deflection  angle  plot  for  the  48-degree  wing/body  AOA, 
Figure  22,  indicates  that  a  canard  deflection  angle  of  -20 
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degrees  could  have  been  chosen,  vice  -25  degrees,  as  the 
canard  deflection  angle  that  nearly  maximized  lift.  This 
choice  would  have  resulted  in  an  canard  absolute  AOA  of  28 
degrees  which  is  much  closer  to  the  average  canard  absolute 
AOA,  for  the  wing/body  angles  of  attack  of  10,  22,  34,  and  40 
degrees,  of  26.8  degrees.  It  appears,  then,  that  the  absolute 
canard  angle  of  attack  to  maximize  lift  enhancement  is 
relatively  constant  and  independent  of  the  wing/body  AOA. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  baseline  study  to  find  the  lift  and  drag  characteristics 
of  a  low-aspect-ratio  wing/body  configuration  across  an  angle- 
of-attack  range  of  -8  to  50  degrees  was  conducted.  A  further 
study  to  find  the  comparative  lift  enhancement  using  the  same 
wing/body  and  a  close-coupled  canard  for  wing/body  angles  of 
attack  of  10,  22,  34,  40,  and  48  degrees  and  canard  deflection 
angles  from  -25  to  25  degrees  was  carried  out.  The  following 
conclusions  were  reached: 

(1)  With  a  canard  properly  located  longitudinally  and 
horizontally,  lift  was  enhanced  at  all  tested  wing/body  angles 
of  attack  when  compared  to  the  baseline  configuration. 

(2)  The  comparative  lift  increase  using  a  canard 
deflection  angle  that  maximized  lift  was  the  most  dramatic  at 
or  near  the  stall  conditions  for  the  wing/body  alone  case.  At 
a  wing/body  AOA  of  22  degrees,  where  the  first  stall  occurred, 
the  lift  was  increased  34%  over  the  baseline  value.  At  a 
wing/body  AOA  of  48  degrees,  where  the  second  stall  occurred, 
the  lift  was  increased  18.3%  over  the  baseline  value.  These 
increases  in  the  lift  are  thought  to  be  due  to  the  vortex  of 
the  canard  constructively  interfering  with  the  main-wing 
vortex,  thereby  delaying  the  onset  of  flow  separation.  Note 
that  these  enhancements  already  take  into  account  the  lift  due 
to  increased  reference  area. 
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(3)  Using  a  canard,  the  lift  was  increased  by  only  3.4% 
for  a  wing/body  AOA  of  10  degrees.  Lacey's  earlier  work 
anticipated  this  result.  Lacey  showed  that  for  angles  of 
attack  less  the  18  degrees  unfavorable  interference  occurred 
between  the  canard  and  wing,  and  little  if  any  lift 
enhancement  could  be  expected  (Ref.  6] . 

(4)  Lift  was  not  enhanced  for  all  canard  deflection 
angles,  which  varied  between  -25  and  25  degrees  in  this  study. 
In  general,  to  maximize  the  lift  the  canard  absolute  angle  of 
attack,  measured  with  respect  to  the  freestream  velocity, 
should  be  about  27  degrees  for  the  tested  configuration.  The 
canard  deflection-angle  range  over  which  lift  was  enhanced  was 
very  much  wing/body  AOA  dependent.  As  the  angle  of  attack 
increased,  the  band  of  canard  deflection  angles  that  lead  to 
lift  enhancement  widened,  but  at  an  AOA  of  10  degrees  this 
band  was  quite  narrow. 

(5)  If  it  is  desired  to  avoid  the  potential  weight  and 
complexity  of  a  moveable  canard  system,  fixing  the  canard  at 
a  -5  degree  deflection  angle  significantly  improved  lift  at 
all  but  the  lowest  wing/body  angles  of  attack.  At  an  AOA  of 
48  degrees,  the  lift  was  only  3%  less  than  the  maximum  lift 
using  a  moveable  canard. 

(6)  Using  a  close-coupled  canard  improved  the  Cj,/Co 
ratio  at  all  angles  of  attack  except  at  10  degrees.  The 
canard  did  cause  a  significant  penalty  of  34.5%  in  the  Cl/Cj, 
ratio  at  an  AOA  of  10  degrees.  It  should  be  noted  that  at 
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this  cruise  flight  condition,  where  the  AOA  of  the  wing  is 
small,  the  AOA  of  the  canard  is  required  to  be  large.  The 
drag  penalty  is  due  to  the  high  induced  drag  of  the  canard. 
At  these  cruise  conditions,  a  high-AOA  lifting  canard  could 
not  be  trimmed,  and  use  of  the  canard  is  a  poor  way  to  achieve 
an  enhanced  lift.  At  this  high-speed  condition,  a  more  likely 
situation  would  be  to  set  the  canard  at  a  neutral  AOA  and  use 
only  the  main  wing  for  lift.  In  this  case,  only  the  main  wing 
would  be  used  as  the  reference  area,  and  a  more  reasonable 
lift  coefficient  could  be  achieved  with  little  drag  penalty. 
In  other  words,  the  canard  is  only  used  at  high  angles  of 
attack,  and  the  apparent  penalty  at  cruise  conditions  is 
provided  by  an  unrealistic  situation. 

Recommendations  for  future  experimental  research  are; 

(1)  Conduct  flow  visualization  tests  of  the  wing/body 
and  canard/wing/body  configuration  to  qualify  and  compare 
flowfields  about  the  model  at  angles  of  attack  of  10,  22,  34, 
40,  and  48  degrees,  and  to  study  the  surface  flow  separation 
and  reattachment  effects.  Perform  wake  surveys  with  a  five- 
hole  probe  at  chosen  angles  of  attack,  with  and  without  the 
canard  at  best  deflection  angle,  to  map  the  location  and 
strength  of  the  vortices. 

(2)  Conduct  a  further  study  using  an  oscillating  canard 
and  the  same  reflection  plane  tunnel  balance,  model  and 
wing/body  angles  of  attack  for  comparison  to  the  static  canard 
configuration. 
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This  thesis  has  shown  that  a  properly  located  close- 
coupled  canard  can  greatly  enhance  lift  at  high  angles  of 
attack  with  no  drag  penalty  when  compared  to  a  wing/body 
configuration.  If  it  is  desired  to  perhaps  more  elegantly 
enhance  lift  at  higher  angles  of  attack  and  avoid  some  of  the 
pitfalls  associated  with  thrust  vectoring  —  such  as  the 
expense,  weight  penalty,  and  excessive  fuel  consumption  —  the 
use  of  a  close-coupled  canard  may  be  an  excellent  choice. 


APPENDIX  A.  MODEL  DESIGN 


The  previous  work  of  Lacey  and  others  was  used  to  set  many 
of  the  design  parameters  of  the  canard/wing  model .  Lacey  used 
a  canard-area-to-wing-area  ratio  of  0.20,  where  the  areas  were 
referenced  to  the  centerline  of  the  fuselage.  It  was 
necessary  that  the  body  of  the  model  used  in  this  design  be 
larger  than  that  used  by  Lacey  to  allow  for  adequate 
instrumentation  space  for  the  canard-moving  mechanism.  Had  a 
similar  area  ratio  to  Lacey's  been  used,  an  exposed  canard 
area  that  was  much  smaller  than  Lacey's  would  have  resulted. 
Accordingly,  an  exposed-canard-area-to-wing  area  (referenced 
to  the  centerline  of  the  fuselage)  ratio  of  0.20  was  chosen, 
which  compared  favorably  to  Lacey's  ratio  of  0.13.  [Ref. 6] 
Aspect  ratios  of  2  for  the  canard  and  3  for  the  wing  were 
used,  based  upon  the  earlier  work  of  Behrbohm.  A  leading  edge 
sweep  of  60  degrees  for  the  canard  and  50  degrees  for  the  main 
wing  were  used  to  ensure  strong  leading  edge  vortices  for  lift 
enhancement .  The  canard  and  wing  were  straight-tapered  and 
taper  ratios  of  .1  and  .15  respectively  were  chosen  based  upon 
existing  aircraft  designs.  Equations  (9),  (10),  and  (11)  were 
used  to  derive  the  dimensions  of  the  planforms.  [Ref.  19] 


AN  »  2 


2? 

C,(l  ♦  k) 


(9) 
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(10) 


MAC  =  -I  ~ 

3  '  C^  *  Ct 

Where : 

AR  Aspect  Ratio 

b  Span 

C,  Length  of  root  chord 

Ct  Length  of  tip  chord 

X  Taper  ratio  C^/C^ 

S  Area  of  wing  or  canard 

MAC  Wing  mean  aerodynamic  chord 

The  airfoil  section  NACA  64A008  was  chosen  for  the  wing 
and  canard  based  upon  Lacey's  previous  work.  Er-El  and 
Seginer  used  flat  plates  for  the  canard  and  wing  of  their 
model  which  resulted  in  sharp  leading  edges  for  the  wing  and 
the  canard  [Ref.  9] .  A  rounded  leading  edge  for  the  wing  and 
canard  was  used  in  this  design  to  more  closely  model  what  is 
found  on  a  number  of  existing  aircraft.  No  attempt  was  made 
to  trip  the  boundary  layer.  The  Reynolds  number  based  on  the 
wing  mean  aerodynamic  chord  was  9.5x10®. 

The  location  of  the  canard  with  respect  to  the  wing  was 
a  critical  dimension.  If  the  canard  was  positioned  too  far 
away  or  too  close,  there  would  be  no  lift  enhancement  or  even 
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possibly  a  lift  degradation.  The  40-percent  exposed  root 
chord  of  the  canard  and  the  quarter-chord  of  the  wing  with 
respect  to  the  centerline  of  the  fuselage  were  the  reference 
points  used  for  the  longitudinal  separation  of  the  canard  and 
wing.  The  ratio  of  the  longitudinal  separation  of  the  canard 
and  wing,  x/c„ac»  should  be  no  greater  than  1.5;  yet,  the 
canard  and  wing  surfaces  should  not  overlap.  Correspondingly, 
an  x/c„ac  =  1-2  was  chosen  which  resulted  in  a  2.33-inch 
separation  between  the  exposed  trailing  edge  of  the  canard  and 
the  exposed  leading  edge  of  the  wing.  Vertically  the  canard 
was  positioned  so  that  the  non-dimensional  distance  of  the 
canard  above  the  wing,  z/c^c  equaled  0.2.  A  value  of  z/c^^  = 
0.2  resulted  in  a  1.9-inch  separation  between  the  canard  and 
wing.  The  pivot  point  of  the  canard  was  40  percent  of  the 
exposed  root  chord.  The  pivot  point  of  the  balance  was  17.18 
inches  from  the  tip  of  the  model.  IRef.  6] 

The  model's  length  of  36  inches,  width  of  4.5  inches, 
height  of  3  inches,  and  semi-span,  measured  from  the 
reflection  plane  to  the  wingtip,  of  12.1  inches  ensured  that 
the  balance  would  be  loaded  by  large  forces  while  the  tunnel 
was  in  operation.  The  size  of  the  model  also  allowed  for 
adequate  space  for  the  canard  positioning  motor.  Initially  it 
was  thought  that  the  canard  would  be  driven  by  remote  radio 
control .  It  was  found  that  the  servo  mechanism  was  unable  to 
hold  the  canard  in  position  at  a  tunnel  AP  of  17  cm  H2O.  A 
lead  screw  mechanism  with  electric  motor  was  then  designed  and 


built  by  lab  personnel  to  positively  control  the  canard  at  all 
tunnel  AP's.  The  electric  cabling  for  the  model  lead  through 
a  hole  in  the  model  base  plate  and  balance  assembly  to  a 
controller  outside  the  tunnel.  The  canard  was  deflected  using 
the  controller  and  a  variable  power  supply.  Lines  drawn  on 
the  canard  section  body  allowed  the  canard  to  be  accurately 
aligned.  Figure  29  gives  the  geometric  characteristics  for 
the  canard  and  wing.  Figures  30  through  36  show  the  model  and 
canard  positioning  motor. 
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Airfoil  Section  (NACA) 

36" 

WING 

64A008 

CANARD 

64A008 

Area  (semi-span),  in'^2 
Projected  Area 

97.9 

50.7 

Exposed  Area 

59.6 

19.3 

Exposed  Semi-Span,  inches 

9.1 

4.4 

Semi-Span,  inches  (centerline) 

12.1 

7.4 

Chord,  inches 

Root  (centerline) 

14 

13.3 

Root  (exposed) 

11 

8 

Tip 

2.1 

0.8 

Aspect  Ratio 

3 

2 

Taper  Ratio 

0.15 

0.1 

Sweepback  Angle,  degrees 
Leading  Edge 

50 

60 

Trailing  Edge 

10.6 

5.5 

Incidence  Angle,  degrees 

0 

Dihedral  Angle,  degrees 

0 

0 

Twist  Angle,  degrees 

0 

- 

Figure  29.  Model  Geometric  Characteristics 
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Figure  33. 
Degrees 


Canard/Wing/Body;  AOA  Approximately  30 
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Figure  35.  Canard  Positioning  Motor  and  Lead-Screw 
Mechanism 
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Figure  36.  Canard/Wing/Body;  Front  View 


APPENDIX  B.  BALANCE  CALIBRATION 


The  external  strain-gage  balance  used  was  built  to 
facilitate  the  measurement  of  normal  and  axial  forces  and 
pitching  moment  in  the  NPS  low  speed  wind  tunnel.  Each 
external  strain-gage  bridge  had  four  active  legs  for  automatic 
temperature  compensation.  The  normal  and  axial  moments  were 
measured  by  two  orthogonal  strain-gage  bridges  cemented  on  the 
balance  column  at  positions  A  and  B  separated  by  a  vertical 
distance  of  26.5  inches,  as  shown  in  Figure  37.  With  the  wind 
tunnel  in  operation  the  force  and  moment  on  the  model  created 
a  different  moment  on  the  upper,  bridge  B,  and  the  lower, 
bridge  A,  strain-gage  bridges.  Once  the  balance  was 
accurately  calibrated,  the  voltage  outputs  from  the  lower  and 
upper  normal  or  axial  bridges  could  be  converted  to  moments 
and  subtracted  from  one  another,  then  divided  by  the  vertical 
separation  of  the  bridges  to  find  the  normal  or  axial  force. 
Figure  38  shows  a  photograph  of  the  balance  and  rotating 
mechanism  used  to  calculate  the  normal  and  axial  forces  and  to 
position  the  model  at  various  angles  of  attack.  Figures  39  to 
41  show  the  wiring  for  the  strain-gage  bridges. 

The  sign  conventions  used  for  the  normal  and  axial  forces 
are  shown  in  Figure  42.  Balance  nomenclature  is  as  follows, 
was  the  voltage  output  from  the  lower  normal  force 
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Balance  Calibration  Rig 
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Figure  38.  Balance  Column  and  Rotating  Mechanism 
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Figure  41.  Strain-Gage  Bridge  Wiring  and  Associated 
Canon  Plugs  that  Connect  to  the  Signal  Conditioners 
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Turntable  at  000  Degrees  Angle  of  Attack 


Normal 

Force 


Axial  Force 


Turntable  at  090  Degrees  Angle  of  Attack 


Normal  Force 


Figure  42.  Normal  and  Axial  Force  Sign  Conventions 
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strain-gage  bridge  and  was  the  voltage  output  from  the 
lower  axial  bridge.  Additionally,  Ebn  was  the  voltage  output 
from  the  upper  normal  bridge  and  was  the  voltage  output 
from  the  upper  axial  bridge.  The  letter  "a"  is  the  vertical 
distance  between  the  upper  bridge  and  the  base  of  the 
calibration  rig  of  3.375  inches.  The  letter  "x"  represents 
the  vertical  distance  between  the  base  of  the  calibration  rig 
and  the  attachment  point  of  the  cable  to  the  rig  of  10.5 
inches.  The  letter  "b"  is  the  vertical  distance  between  the 
upper  and  lower  strain-gage  bridges  of  26.5  inches. 

Equations  (12)  and  (13)  are  the  basic  equations  used  to 
find  the  normal  and  axial  moments  on  the  balance  at  the  lower 
and  upper  bridge  positions.  For  equation  (12),  the  electri¬ 


An 

*12 

Mb. 

(12) 

■*21 

*22 

Bii 

5i2 

Mbs 

(13) 

Ebb 

^21 

B22 

Mbb 

cal  output  in  volts  at  the  lower  normal  and  axial  strain-gage 
bridges,  bridge  location  A,  equalled  a  coefficient  matrix.  A, 
multiplied  by  a  moment  matrix  in  ft-lbs  for  the  lower  bridge 
that  consisted  of  normal  and  axial  components.  For  equation 
(13),  the  electrical  output  at  the  upper  normal  and  axial 
strain-gage  bridges,  bridge  location  B,  equalled  a  coefficient 
matrix,  B,  multiplied  by  a  moment  matrix  for  the  upper  bridge 
that-  consisted  of  a  normal  and  axial  component.  The 
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coefficient  matrices  were  found  by  calibrating  the  balance 
with  known  weights.  In  general,  the  diagonal  components  of 
the  A  and  B  coefficient  matrices,  A^^  Ki2i  Bn,  and  832,  were 
larger  than  the  off-diagonal  components  by  at  least  one  or  two 
orders  of  magnitude.  The  off-diagonal  components  were  due  to 
the  slight  misalignment  of  the  strain  gages  which  caused 
interaction  between  them.  Under  ideal  conditions,  the  off- 
diagonal  components  are  close  to  zero.  In  general,  balance 
interactions  are  due  to  misalignment  of  the  balance  or  plastic 
deformation  of  the  balance  [Ref.  18] .  Balance  misalignment 
causes  a  linear  or  first  degree  error  whereas  balance  plastic 
deformation  would  cause  a  second  order  or  non-linear  error 
[Ref.  18] .  The  calibration  procedure  used  only  corrected  for 
misalignment  of  the  strain  gages,  for  the  balance  was  never 
plastically  deformed.  Equations  (14)  through  (17)  were  used 
during  the  calibration  to  calculate  axial  and  normal  moments 
at  the  upper  and  lower  bridges,  using  the  moment  arms  a,  b, 
and  X. 

-  Normal  Force  (b  +  a  +  x)  (14) 

=  Axial  Force  (b  +  a  +  x)  (15) 


Mgf,  =  Normal  Force  {a  +  x) 


(16) 


Mgj^  =  Axial  Force  la  +  x) 


(17) 


Before  the  balance  could  be  calibrated  a  loading  rig  had 
to  be  designed.  One  of  the  main  requirements  was  that  it  be 
very  rigid  to  preclude  any  chance  of  elastic  or  plastic 
deformation  of  the  rig.  Secondly,  it  had  to  have  arms  that 
could  rotate  and  move  vertically  for  attaching  the  cabling  and 
weights  so  that  the  balance  could  be  calibrated  in  various 
directions  with  different  moment  arms.  The  calibration  rig  was 
pulled  horizontally  using  known  weights  suspended  on  a  cable, 
stand,  and  pulley  assembly.  The  calibration  rig  is  shown  in 
Figures  43  to  45. 

Prior  to  calibrating  the  balance  the  span  voltages  for 
the  bridges  were  set .  They  were  set  by  suspending  a  weight  of 
about  10  lbs  from  the  calibration  rig  and  turning  the  span 
voltage  adjustment  knob  for  each  signal  conditioi»er  so  that 
the  readout  on  a  Hewlett-Packard  digital  multimeter  was  about 
1  volt.  Once  the  span  voltages  were  set,  they  were  not 
altered  for  the  duration  of  the  experiment. 

The  diagonal  components  and  Bu  were  found  by  rotating 
the  balance  to  000  degrees  angle  of  attack  for  loading  the 
balance  in  the  pure  normal  direction,  then  loading 
successively  larger  weights  in  approximately  2-lb  increments 
on  the  calibration  rig  and  plotting  the  resultant  voltage 
output  versus  the  weight.  Figures  46  and  47  and  Table  4  show 
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plots  of  the  magnitude  of  voltage  output  in  volts  versus 
weight  in  pounds  and  the  data  used  for  the  plots.  The 
gradients  of  the  voltage  versus  weight  plots,  E^j/N  and  Ebn/N, 
from  Figures  46  and  47  respectively,  were  used  in  equations 
(18)  and  (19)  to  find  the  diagonal  components  and  Bn  of 
the  coefficient  matrices. 


E^/ Normal  Force 
”  (h+a+x) 


(18) 


E^/ Normal  Force 
“  "  - Ta+x) - 


(19) 


The  diagonal  components  A22  and  B22  were  found  by  rotating 
the  balance  to  90  degrees  AOA  for  loading  the  balance  in  the 
pure  axial  direction,  then  following  a  similar  procedure  to 
that  used  to  find  An  and  Bn.  The  gradients  used  for  voltage 
output  versus  weight,  E;^;^/A  and  E^^/A,  are  found  in  Figures  48 
and  49.  Table  5  gives  the  data  used  in  the  plots.  The 
gradients  were  used  in  equations  (20)  and  (21)  to  find  the 
diagonal  components  A22  and  B22. 


Ejj^/ Axial  Force 
(b+a+x) 


(20) 
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E  B  AXIAL  (V) 


Figure  49.  vs.  Weight 


94 


Table  4.  Strain-Gage  Bridge  Outputs  for  Pure  Normal 
Forces 


EA  AXIAL  (V) 

EB  AXIAL  (V) 

EA  NORM  (V) 

EB  NORM  (V) 

1.0754 

-0.005 

-0.026 

0.146 

0.097 

2.0714 

-0.006 

-0.027 

0.266 

0.199 

3.0780 

-0.008 

-0.029 

0.387 

0.299 

4.2155 

-0.010 

-0.032 

0.535 

0.421 

5.2144 

-0.012 

-0.034 

0.660 

0.520 

6.2190 

-0.013 

-0.037 

0.782 

0.630 

8.2180 

-0.017 

-0.044 

1.013 

0.817 

10.2215 

-0.019 

-0.051 

1.261 

1.025 

13.2216 

-0.024 

-0.063 

1.675 

1.362 

17.2241 

-0.032 

-0.082 

2.139 

1.750 

21.2328 

-0.038 

-0.099 

2.607 

2.139 

0.0000 

-0.002 

-0.028 

0.011 

-0.012 

Table  5.  Strain-Gage  Bridge  Outputs  for  Pure  Axial 
Forces 


HKSIil 

EA  AXIAL  (V) 

EB  AXIAL  (V) 

EB  NORM  (V) 

1.0754 

0.109 

0.092 

0.002 

-0.004 

2.0714 

0.215 

0.161 

0.004 

-0.002 

3.0780 

0.317 

0.229 

0.005 

-0.002 

4.2155 

0.441 

0.309 

0.005 

-0.002 

5.2144 

0.549 

0.380 

0.007 

-0.002 

6.2190 

0.635 

0.435 

0.008 

-0.001 

8.2180 

0.845 

0.577 

0.009 

-0.001 

10.2215 

1.039 

0.703 

0.013 

0.000 

13.2216 

1.352 

0.910 

0.017 

0.000 

17.2241 

1.770 

1.181 

0.022 

0.001 

21.2328 

2.181 

1.453 

0.025 

0.003 

0.0000 

-0.001 

0.016 

0.000 

-0.004 
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(21) 


(Egj^/ Axial  Force) 

B,2  =  -  ■  - 

“  a+x 

The  off-diagonal  coefficients  A21/  B^j/  and  B21  were 

found  by  choosing  an  angle  of  attack  and  a  weight  for  a 
combined  loading, i.e,  a  turntable  angle  of  attack  equal  to  55 
degrees  and  a  weight  of  17  pounds,  then  using  the  resultant 
electrical  output  and  known  diagonal  coefficients  to  solve  the 
equation  for  the  off-diagonal  coefficient.  Equation  (22)  is 
an  example  of  solving  for  an  unknown  off-diagonal  coefficient. 


^12  ~ 


-  All  (Normal  Force)  (a+b+x) 


(Axia. 


^orce)  la+b+x) 


Once  all  the  calibration  coefficients  were  found  for  the 
matrices,  the  A  or  B  matrix  was  then  inverted  and  an  unknown 
moment  out  of  the  normal  or  axial  bridges  could  be  solved  for. 
Equations  (23)  and  (24)  are  the  final  calibration  matrix 
equations  for  the  upper  and  lower  strain-gage  bridges.  To 
find  the  normal  force,  the  moment  from  the  upper  normal  bridge 
at  position  B  was  subtracted  from  the  moment  from  the  lower 
normal  bridge  at  position  A.  The  result  was  then  divided  by 
b,  the  distance  between  the  strain-gage  bridges.  This  method 
subtracted  out  any  residual  moment  on  the  model  for  the 
determination  of  moments  due  solely  to  normal  and  axial 
forces.  The  axial  force  was  found  in  a  similar  fashion.  See 
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equations  (25)  and  (26)  .  The  lift  and  drag  were  computed 
using  equations  (3)  and  (4)  found  in  the  Apparatus  section. 


I 


"an 

27.26281 

-.48797 

Waa 

-.04932 

32.53814 

11.30306  -.70957 

Ebn 

«BA 

-.1350659  17.25848 

Ea. 

(23) 


(24) 


Normal  Force  = 


- 


Wo, 


(25) 


Axial  Force  = 


Waa  -  «aA 


(26) 


C 
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APPENDIX  C.  DATA  ACQUISITION  PROGRAM 


0100  'THIS  DATA  ACQUISTION  PROGRAM  HAS  BEEN  MODIFIED  BY  LT . 
YUAN  AND  LT  KERSH 

0110  'FOR  THE  EXTERNAL  STRAIN  GAGE  BALANCE.  THE  USER  IS 
ENCOURAGED  TO  MODIFY 

0120  'THE  ORIGINAL  PROGRAM,  READ. BAS  FOR  THEIR  OWN  PURPOSE. 
1000  DEF  SEG:  CLEAR,  &HFE00:  GOTO  1030  'Begin  PCIB 

Program  Shell 

1010  GOTO  2900  'User  program 

1020  GOTO  2670  'Error  handling 

1030  I=&HFE00  'Copyright  Hewlett-Packard 

1984,1985 

1040  PCIB. DIR$=ENVIRON$ ("PCIB") 

1050  I$*PCIB.DIR$+"\PCIBILC.BLD" 

1060  BLOAD  I$,I 

1070  CALL  I (PCIB.DIRS, 1%, J%) :  PCIB.SEG=I% 

1080  IF  J%=0  THEN  GOTO  1120 
1090  PRINT  "Unable  to  load."; 

1100  PRINT  "  (Error  #";J%;")" 

1110  END 
1120  ' 

1130  DEF  SEG=PCIB.SEG:  O.S=5:  C.S=10;  I.V»15 

1140  I.C»20:  L.P=25:  LD.FILE=30 

1150  GET.MEM=35:  L.S=40:  PANELS=45:  DEF.ERR=50 

1160  PCIB.ERR$=STRING$(64,32) ;  PCIB.NAME$=STRING$ (16, 32) 

1170  CALL  DEF. ERR (PCIB. ERR, PCIB. ERR$, PCIB. NAMES, PCIB. GLBERR)  : 

PCIB.BASERR*255 

1180  ON  ERROR  GOTO  1020 

1190  J=-l 

1200  I$=PCIB.DIR$+"\PCIB.SYN" 

1210  CALL  O.S(I$) 

1220  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1230  1=0 

1240  CALL 

I .V (I, READ . REGISTER, READ . SELFID, DEFINE, INITIALIZE . SYSTEM) 
1250  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1260  CALL 

I .V ( I, ENABLE . SYSTEM, DISABLE . SYSTEM, INITIALIZE, POWER. ON) 

1270  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1280  CALL  I. V (I, MEASURE, OUTPUT, START, HALT) 

1290  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1  3  0  0  CAL 

I .V ( I, ENABLE . INT . TRIGGER, DISABLE . INT . TRIGGER, ENABLE . OUTPUT, 

I SABLE. OUTPUT) 

1310  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1320  CALL  I. V(I, CHECK. DONE, GET. STATUS, SET. FUNCTION,SET. RANGE) 
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1330  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1340  CALL  I. V(I, SET. MODE, WRITE. CAL,READ. CAL, STORE. CAL) 

1350  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1360  CALL  I. V( I, DELAY, SAVE. SYSTEM, J,J) 

1370  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1380  1=1 

1390  CALL  I.V(I, SET. GATETIME,SET. SAMPLES, SET. SLOPE, SET. SOURCE) 

1400  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1410  CALL  I. C( I, FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER) 

1420  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1430  CALL  I. C (I, INTERVAL, RATIO, TOTALIZE, RIOOMILLI) 

1440  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1450  CALL  I.C(I,R1,R10,R100,R1KILO) 

14  60  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1470  CALL  I. C( I, R1 OMEGA, RIO OMEGA, CHAN. A, CHAN. B) 

1480  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1490  CALL  I. C( I, POSITIVE, NEGATIVE, COMN, SEPARATE) 

1500  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1510  1=2 
1520  1=3 

1530  CALL  I. V(I, ZERO. OHMS, SET. SPEED, J,J) 

1540  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1550  CALL  I.C(I,DCVOLTS,ACVOLTS,OHMS,R200MILLI) 

1560  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1570  CALL  I.C(I,R2,R20,R200,R2KILO) 

1580  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1590  CALL  I.C(I,R20KILO,R200KILO,R2MEGA,R20MEGA) 

1600  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1610  CALL  I.C(I,AUTOM,R2.5,R12.5, J) 

1620  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1630  1=4 

1640  CALL 

I .  V  { I ,  SET . COMPLEMENT, SET . DRIVER, OUTPUT . NO . WAIT, ENABLE . HANDSH 
AKE) 

1650  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1660  CALL 

I . V ( I , D I SABLE . HANDSHAKE , SET . THRESHOLD , SE T . S TART .BIT,SET.NUM. 
BITS) 

1670  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1680  CALL  I.V(I,SET.LOGIC.SENSE, J, J, J) 

1690  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1700  CALL  I. C  (I, POSITIVE, NEGATIVE, TWOS, UNSIGNED) 

1710  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1720  CALL  I.C(I,OC,TTL,RO,R1) 

1730  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1740  CALL  I.C(I,R2,R3,R4,R5) 

1750  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1760  CALL  I.C(I,R6,R7,R8,R9) 

1770  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1780  CALL  I.C(I,R10,R11,R12,R13) 

1790  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
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1800  CALL  I.C(I,R14,R15,R16,  J) 

1810  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1820  1=6 

1830  CALL 

I .  V  ( I , SET . FREQUENCY , SET . AMPLITUDE, SET . OFFSET , SET . SYMMETRY ) 
1840  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1850  CALL  I.V(I,  SET.  BURST.  COUNT,  J,J,J) 

1860  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1870  CALL  I. C( I, SINE, SQUARE, TRIANGLE, CONTINUOUS) 

1880  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1890  C;a.L  I.C(I,GATED,BURST,  J,  J) 

1900  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1910  1=7 

1920  CALL 

I  .V  (I, AUTOSCALE, CALIBRATE, SET . SENSITIVITY, SET .VERT . OFFSET) 
1930  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1940  CAL  L 

I  .V  (I,  SET .COUPLING, SET .POLARITY, SET . SWEEPSPEED, SET .DELAY) 
1950  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1960  CALL 

1  .V(I,  SET. TRIG. SOURCE, SET. TRIG. SLOPE, SET . TRIG . LEVEL, SET. TRIG 
.MODE) 

1970  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

1980  CALL 

I ,  V  ( I ,  GET . SINGLE .WF, GET . TWO .WF, GET .VERT . INFO, GET . TIMEBASE . IN 
FO) 

1990  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2000  CALL 

I .  V  ( I ,  GET . TRIG . INFO, CALC . WFVOLT, CALC . WFTIME, CALC . WF . STATS ) 
2010  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2020  CALL 

I .  V  ( I ,  CALC . RI SETIME , CALC . FALLTIME , CALC . PERI OD , CALC . FREQUENCY 

) 

2030  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2040  CALL 

I .  V  ( I ,  CALC . PLUSWIDTH, CALC . MINUSWIDTH, CALC . OVERSHOOT , CALC . PRE 
SHOOT) 

2050  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2  0  6  0  CAL 

I . V  ( I ,  CALC . PK . TO . PK, SET . TIMEOUT, SCOPE . START , MEASURE . SINGLE . 
F) 

2070  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2080  CALL  I.V(I,MEASURE.TWO.WF, J, J, J) 

2090  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2100  CALL  I.C(I,R10NANO,R100NANO,R1MICRO,R10MICRO) 

2110  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2120  CALL  I.C(I,R100MICRO,R1MILLI,R10MILLI,R100MILLI) 

2130  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2140  CALL  I.C(I,R1,R10,R20NANO,R200NANO) 

2150  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2160  CALL  I.C(I,R2MICRO,R20MICRO,R200MICRO,R2MILLI) 


100 


S  f 


2170  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2180  CALL  I.C(I,R20MILLI,R200MILLI,R2,R20) 

2190  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2200  CALL  I.C(I,R50NANO,R500NANO,R5MICRO,R50MICRO) 

2210  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2220  CALL  I .C (I, R500MICRO, R5MILLI, R50MILLI, R500MILLI ) 

2230  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2240  CALL  I .C (I,R5,R50,CHAN.A,CHAN.B) 

2250  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2260  CALL  I .C (I, EXTERNAL, POSITIVE,NEGATIVE, AC) 

2270  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2280  CALL  I .C (I, DC, TRIGGERED, AUTO.TRIG,AUTO.LEVEL) 
2290  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2300  CALL  I. C( I, XI, XI 0, STANDARD, AVERAGE) 

2310  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2320  1=8 

2330  CALL  I .V(I,OPEN.CHANNEL,CLOSE.CHANNEL, J, J) 

2340  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2350  CALL  C.S 

2360  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2370  I$=PCIB.DIR$+"\PCIB.PLD" 

2380  CALL  L.P(I$) 

2390  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2400  I$="DMM.01":  1=3:  J=0:  K=0:  L=1 
2410  CALL  DEFINE(DMM.01,I$,I, J,K,L) 

2420  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2430  I$="Func.Gen.01":  1=6:  J=0:  K=l:  L=1 
2440  CALL  DEFINE (FUNC .GEN. 01, 1$, I, J, K,  L) 

2450  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2460  I$="Scope.01":  1=7:  J=0:  K=2:  L=1 
2470  CALL  DEFINE (SCOPE . 01, 1$, I, J, K, L) 

2480  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2490  I$="Counter.01":  1=1:  J=0:  K=3:  L=1 
2500  CALL  DEFINE (COUNTER. 01, 1 $, I,  J,K,L) 

2510  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2520  I$="Dig.In.01":  1=4:  J=0:  K=4:  L=1 
2530  CALL  DEFINE (DIG. IN. 01, 1$, I, J,K,L) 

2540  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2550  I$="Dig.Out.01":  1=4:  J=l:  K=4:  L=1 
2560  CALL  DEFINE (DIG.OUT.Ol, 1$, I, J,K,L) 

2570  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2580  I$="Relay.Act.0l":  1=8:  J=0:  K=5:  L=1 
2590  CyOiL  DEFINE  (RELAy.ACT.01,I$,  I,  J,K,L) 

2600  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2610  I$="Relay.Mux.01":  1=2:  J=0:  K=6:  L=1 
2620  CALL  DEFINE (RELAY. MUX. 01, 1$, I, J,K, L) 

2630  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
2640  I$=ENVIRON$ ("PANELS") +" \PANELS.EXE" 

2650  CALL  L.S(I$) 

2660  GOTO  1010 

2670  IF  ERR=PCIB.BASERR  THEN  GOTO  2700 
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2680  PRINT  "BASIC  error  #";ERR;"  occurred  in  line  ";ERL 
2690  STOP 

2700  TMPERR=PCIB.ERR:  IF  TMPEKR=0  THEN  TMPERR=PCIB.GLBERR 
2710  PRINT  "PC  Instrument  error  #";TMPERR;"  detected  at  line 
";ERL 

2720  PRINT  "Error:  ";PCIB.ERR$ 

2730  IF  LEFTS (PCIB.NAME$,1)<>CHR$ (32)  THEN  PRINT  "Instrument: 
";PCIB. NAMES 
2740  STOP 

2750  COMMON  PCIB.DIRS/PCIB.SEG 

2760  COMMON  LD.FILE, GET. MEM, PANELS, DEF. ERR 

2770  COMMON 

PCIB . BASERR, PCIB . ERR, PCIB . ERRS , PCIB . NAMES , PCIB . GLBERR 
2780  COMMON 

READ .REGISTER, READ . SELF ID, DEFINE, INITIALIZE . SYSTEM, ENABLE . SY 
STEM, DISABLE. SYSTEM, INITIALIZE,POWER. ON, MEASURE, OUTPUT, START 
, HALT, ENABLE . INT . TRIGGER, DISABLE . INT . TRIGGER, ENABLE . OUTPUT , D 
I SABLE . OUTPUT, CHECK . DONE, GET . STATUS 

2790  COMMON 

SET . FUNCTION, SET . RANGE, SET .MODE, WRITE . CAL, READ . CAL, STORE . CAL 
, DELAY, SAVE . SYSTEM, SET .GATETIME, SET . SAMPLES, SET . SLOPE, SET . SO 
URGE, ZERO . OHMS , SET . SPEED, SET . COMPLEMENT, SET . DRIVER, OUTPUT . NO 
. WAI T , ENABLE . HANDSHAKE , D I SABLE . HANDSHAKE 

2800  COMMON 

SET . THRESHOLD , SET . START . BIT , SET . NUM . BITS , SET . LOGIC . SENSE , SET 
.FREQUENCY, SET .AMPLITUDE, SET .OFFSET, SET . SYMMETRY, SET .BURST . C 
OUNT,  AUTOSC/ILE,  CALIBRATE,  SET .  SENSITIVITY,  SET  .VERT .  OFFSET,  SET 
. COUPLING, SET .POLARITY, SET . SWEEPSPEED 

2810  COMMON 

SET . DELAY, SET . TRIG . SOURCE, SET . TRIG . SLOPE, SET . TRIG . LEVEL, SET . 
TRI G . MODE , GET . S INGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . T IMEB?  S 
E . INFO, GET . TRIG . INFO, CALC . WFVOLT, CALC . WFTIME, CALC . WF . STATS, C 
ALC . RI SET IME , CALC . FALLTIME , CALC . PERI OD 

2820  COMMON 

CALC .FREQUENCY, CALC .PLUSWIDTH, CALC.MINUSWIDTH, CALC . OVERSHOOI 
, CALC . PRESHOOT, CALC . PK . TO . PK, SET . TIMEOUT, SCOPE . START, MEASURE 
. S INGLE . WF , MEASURE . TWO . WF , OPEN . CHANNEL , CLOSE . CHANNEL 
2830  CO  MMON 

FREQUENCY, AUTO .FREQ, PERIOD, AUTO .PER, INTERVAL, RATIO, TOTALIZE, 
Rl 0  OMILLI , Rl , Rl 0 , Rl 0  0 , RIKILO, Rl OMEGA, Rl 0  OMEGA, CHAN .A, CHAN . B , 
POSITIVE, NEGAT I VE , COMN , SEPARATE , DCVOLTS , ACVOLT S , OHMS , R2  0  OM I L 
LI ,  R2 , R2 0 , R2 0 0 , R2KILO, R2 OKILO, R2 0 OKILO 

2840  COMMON 

R2MEGA,  R20MEGA,  AUTOM,  R2 . 5,  R12 . 5, POSITIVE,  NEGATIVE,  TWOS,  UNSIG 
NED,OC,TTL,R0,Rl,R2,R3,R4,R5,R6,R7,R8,R9,R10,Rll,R12,R13,R14 
, Rl 5 , Rl 6 , S INE , SQUARE , TRIANGLE , CONTINUOUS , GATED , BURST , Rl ONANO 
, RIO ONANO, RIMICRO, RIOMICRO, RIOOMICRO 

2850  COMMON 

RIMILLI ,  Rl OMILLI, RIO OMILLI , Rl, RIO, R2 ONANO, R200NANO, R2MICRO, R 
20MICRO,R200MICRO,R2MILLI,R20MILLI,R200MILLI,R2,R20,R50NANO, 
R500NANO,R5MICRO,R50MICRO,R500MICRO,R5MILLI,R50MILLI,R500MIL 
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LI , R5 , R50 , CHAN .A, CHAN . B, EXTERNAL, POSITIVE 

2860  COMMON 

NEGATIVE, AC, DC, TRIGGERED, AUTO . TRIG, AUTO . LEVEL, XI , XI 0 , STANDAR 
D, AVERAGE 

2870  COMMON 

DMM. 01, FUNC. GEN. 01, SCOPE. 01, COUNTER. 01, DIG. IN. 01, DIG. OUT. 01, 
RELAY . ACT . 0 1 , RELAY . MUX . 0 1 
2880  'End  PCIB  Program  Shell 
2890  ' 

2900  'Program  to  scan  with  the  DMM  and  RELAY. MUX. 01 

2910  'This  program  was  writen  by  T.SESTAK  and  modified  by 

2920  'P.  ROANE  ,  P.  RABANG,  J. SOMMERS,  and  J.  Kersh  for  use 

2925  'with  the  balance  designed  by  Prof.  Schmidt  et  al. 

2926  ' 

2930  ' 

2940  'The  next  section  after  the  SHELL  program  directs  reading 
2950  'the  voltages  from  the  balance,  computes  forces  measured 
2960  'by  the  strain  gages,  then  stores  the  values  in  two 
arrays, 

2970  'one  for  the  TARE  one  for  FORCE.  This  data  file  can  then 
2980  'be  used  for  graphs  or  other  displays.  Each  test  run 
2990  'will  create  a  file  on  the  A  and  C  drives  that  consist 
3000  'raw  data  that  has  been  manipulated  into  normal,  axial, 
3010  'lift,  and  drag  forces.  The  raw  data  is  stored  on  the 
3020  'C  drive  in  a  BALANFILE. 

3030  'Arrays  are  dimensioned  here 

3  0  4  0  DIM 

READ ING t 7 ] , FORCE t 1 4 0 , 9 ] , TARE [ 8 ] , TREAD t 7 , 1 0 ] , LAB t 7 ] , DROP [ 7 ] 
3050  COLOR  14,1,1 
3060  CLS 
3070  KEY  OFF 

3080  AOA=0  'Last  modified  on  2 7 Sep 90 

3085  TEMP=0 
3090  VALUE=5 

3092  CLSrLOCATE  11,10:PRINT  "IN  THIS  PROGRAM," 

3094  LOCATE  12,10:PRINT  "YOU  HAVE  TO  ANSWER  ALL  THE  QUESTION 
BY  CAPITAL  LETTERS." 

3096  LOCATE  13,10:PRINT  "SO,  PLEASE  TURN  ON  THE  'CAPS  LOCK', 
THANKS ! " 

3098  LOCATE  15, 10: INPUT  "ENTER  <CR>  TO  CONTINUE"; INPT$ 

3100  CNDAOA=0 

3110  CLSzLOCATE  11, 28 :PRINT"SETTING  UP  DATA  FILES" 

3115  LOCATE  13, 20: INPUT  "ENTER  THE  OUTPUT  DATA  FILE  NAME";D$ 
3120  D$=D$+".DAT" 

3130  'The  program  writes  the  data  to  several  files. 

3140  STATEFILE$  «  "C:\PCIB\WIND.HPC"  'stored  in  PCIB 

subdirectory 

3150  DATAFILES  =  "C: \KERSH\"+D$  'stored  on  drive 

C 

3160  DISKFILES  *  "A:"+D$  'stored  on  drive 

A 
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3170  BALANFILE$  *=  "C :  \KERSH\B"+D$ 

C 

3180  ' 

3190  RELAY. SETTLING. TIME  =  .90  'Can  be  changed,  but 

don't ! 

3200  LOCATE  16, 35  :PRINT’’D  ONE" 

3210  CALL  DELAY (VALUE) 

3220  ' 

3230  CLSrLOCATE  12,28 :PRINT"INITIALIZING  INSTRUMENTS" 

3240  CALL  INITIALIZE. SYSTEM (STATEFILE$) 

3250  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
3260  CALL  ENABLE . SYSTEM 

3270  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
3280  LOCATE  16, 35 :PRINT"D  ONE" 

3290  CALL  DELAY (VALUE) 

3300  ' 

3310  'This  part  of  the  program  is  to  preserve  the  data  if 
3320  'if  the  program  is  aborted  in  mid-run.  Parity  errors 
3330  ' in  the  Hewlett  Packard  PC  Instruments  setup  caused  by 
3340  'electrical  noise  and  undervoltage  at  NPS  requires 
3350  'this.  A  voltage  regulated,  uninteruptible  power  supply 
3360  'would  ameliorate  this  problem.  Just  in  case-  this  little 
3370  'sequence  allows  reentry  into  the  program  and  the  data 
3380  'arrays  with  minimal  inconvenience. 

3390  ' 

3400  CLS:LOCATE  12, 20 ; INPUT"WERE  YOU  INTERRUPTED  (Y  OR  N)";A$ 
3410  IF  A$="y"  THEN  GOTO  3500 
3420  ' 

3430  'The  next  two  variables  are  counters  in  the  arrays 
3440  'FORCE  and  TARE 
3450  ' 

3460  TRIAL  »  0 
3470  TRY  =  0 
3480  GOTO  3690 
3490  ' 

3500  LOCATE  14, 15; INPUT  "WHAT'S  THE  INTERRUPTED  FILE 
NAME";ITDF$ 

3510  ITDF$«"C;\LAWRENCE\"+ITDF$+".DAT" 

3520  OPEN  ITDF$  FOR  INPUT  AS  #1 

3530INPUT#1 
TARE  (1) ,  TARE  (2) ,  TARE  (3) ,  TARE  (4) ,  TARE  (5) ,  TARE  (6) ,  TARE  (7) ,  TARE 
(8) 

3540  FOR  X  *=  1  TO  140 

3550  INPUT  #1, 

FORCE (X, 1) , FORCE (X, 2) , FORCE (X, 3) , FORCE (X, 4) , FORCE (X, 5) , FORCE 
(X, 6) , FORCE (X, 7) , FORCE (X, 8) , FORCE (X, 9) 

3560  IF  FORCE (X,1)«0  THEN  ABCD»X;GOTO  3590 

3570  AOA»FORCE(X,2) 

3575  TEMP-FORCE (X, 9) 

3580  NEXT  X 
3590  TRY=ABCD-1 
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3600  TRIAL*ABCD-1 
3610  CLOSE  #1 
3620  GOTO  3700 
3640  ' 

3650  'A$  is  used  as  a  marker  for  interrupted  run  sequences 
3660  'in  the  program,  it  is  set  to  "N"  so  the 
3670  'uninterrupted  sequences  are  used  unless  otherwise 
directed 


3680  ' 

3690  A$=’’N« 

3700  ' 

3710  ' 

3720  'Prompt  to  begin  each  scan  or  quit  program  if  desired 
3730  ' 

3740  CLS:LOCATE  12,10 

3750  INPUT  "TO  START  SCAN  ENTER  ANY  KEY  EXCEPT  Q,  Q  TO 
QUIT" /ANSWERS 

3760  IF  ANSWERS  =  "Q"  THEN  GOTO  6655 
3770  ' 

3780  'This  section  enters  AOA,  Canard  AOA,  and  Temperature  for 
each 

3790  'trial  and  displays  them  in  the  printout. 

3800  CLS: LOCATE  12,10 

3810  PRINT  "THE  LAST  ANGLE  OF  ATTACK  WAS  ";AOA 
3820  LOCATE  13, 10; PRINT  "THE  LAST  CANARD  AOA  WAS  "/CNDAOA 
3830  LOCATE  14, 10 .-PRINT  "THE  LAST  TEMPERATURE  WAS  "/TEMP 
3840  LOCATE  16, 10: INPUT  "ENTER  THE  ANGLE  OF  ATTACK  (AOA)  FOR 
THIS  TRIAL" /AOA 

3850  LOCATE  17, 10;  INPUT  "ENTER  THE  CANARD  AOA  FOR  THIS 
TRIAL" /CNDAOA 

3860  LOCATE  18, 10: INPUT  "ENTER  THE  TEMPERATURE  FOR  THIS 

TRIAL" /TEMP 

3870  READING (l)*AOA 

3880  AOA»  90  -  AOA  'Model  mounted  at  a  turntable  AOA  of 


090  deg 

3890  'OAA*  'Voltage  offset  values  can  be  entered 

here 

3900  'OBA=  'if  the  offset  values  are  not  time  or 

3910  'OAN«  'temperature  dependent. 

3920  'OBN= 

3930  CYCLE  =  0 

3940  CLS;LOCATE  13,10;INPUT  "THE  Eo  FOR  EAA  IS"/OAA 

3950  CLS: LOCATE  14, 10: INPUT  "THE  Eo  FOR  EBA  IS"/OBA 

3960  LOCU^TE  15, 10: INPUT  "THE  Eo  FOR  EAN  IS"/OAN 

3970  LOCATE  16,10:INPUT  "THE  Eo  FOR  EBN  IS"/OBN 

3980  'The  user  may  reject  the  current  readings  and  input  a  new 

set. 


3990  ' 

4000  CLS 

4010  PRINT"********************  DIRECT  BALANCE  READINGS 

*********************** 
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4020  PRINT"  CHECK  OF  SYSTEM  OPERATION 

4030  PRINT 

4031  PRINT"IN  VOLTS 

4032  ' 

4040  PRINT"******** 

4041  ' 

4060  'This  file  is  for  storing  the  direct  voltage  readings  and 
averages . 

4070  'The  data  file  is  continually  appended  for  each  tunnel 
run. 

4080  'The  data  is  for  further  analysis  of  the  direct  voltage 
readings . 

4090  OPEN  BIVLANFILES  FOR  APPEND  AS  #3 
4100  ' 

4110  FOR  CNT  =  1  TO  10 

4120  FOR  CHANNEL  =  2  TO  5 

4130  CALL  OUTPUT (RELAY. MUX. 01,  CHANNEL) 

4140  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

4150  CALL  DELAY (RELAY. SETTLING. TIME) 

4160  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

4170  CALL  MEASURE (DMM. 01,  READING [CHANNEL] ) 

4180  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

4185  READING (CHANNEL) ^READING (CHANNEL) 

4190  TREAD (CHANNEL, CNT)  «  READING (CHANNEL) 

4200  NEXT  CHANNEL 
4210  PRINT  USING 

"  +.««««««  +.««««««  +.«»»»«» 
+  .######  +.###### 

+  . ######";READING (2) , READING (3) , READING (4) , READING  (5) 

4220  PRINT  #3,  USING 

"  +###.#  +.######  +.«#«#«#  +.«««#«« 
+  .######  +.###### 

-•-.######";READING(1)  ,READING(2)  ,READING(3)  ,READING(4)  ,R£ADIN 
G(5) 

4230  NEXT  CNT 
4240  ' 

4250  '  CALL  SUBROUTINE  TO  AVERAGE  READINGS 
4260  GOSUB  6690 
4270  ' 

4280  PRINT" - 


4290  PRINT  USING 

"MEAN  VALUE  4. ###««#  4. «####«  +.«««««« 

+  .««««««  +.###### 

-t- .  ######";READING  (2) ,  READING  (3)  ,  READING  (4) ,  READING  (5) ,  READIN 
G(  6),  READING  (7) 

4300  PRINT  «3,  USING 

"  +###.#  4. ««««««  +.««««««  +.«««««« 

+  .«««##«  +.###### 

+ .######";READING (1) , READING (2) , READING (3) , READING (4) , READIN 
G(5) , READING ( 6) , READING (7) 


EAA 

*** 


EBA 


•kltlf 


EAN 
*  *  * 


EBN" 

***•• 
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4310  CLOSE  #3 
4320  PRINT"  ".‘BEEP 

4330  PRINT"<CR>  TO  CONTINUE,  "1"  TO  GET  NEW  READINGS" 

4340  INPUT  XYZ 

4350  IF  XyZ=l  GOTO  3940 

4360  ' 

4370  'The  voltage  readings  from  the  balance  are  corrected  with 
the  zero  offset 

4380  'so  that  their  slope  passes  through  the  coordinate 

origin. 

4390  'The  corrected  values  are  then  applied  to  the  balance 

interaction 

4400  'equations. 

4410  ***************  CONVERT  SIGNAL  TO  FORCES  *************** 
4420  r*************************************^-^  **************** 
4430  'The  calibration  matrices  were  multiplied  by  a  factor  of 
1000. 

4440  'Ensure  that  a  gain  of  1000  is  set  on  the  Pacific 

Amplifier. 

4450  ' 

4460  ' 

4470  'The  next  section  was  added  to  compensate  for  the 

dependence 

4480  'of  Channel  3,  strain  bridge  EBA,  on  temperature. 

Channel 

4490  '3  was  found  to  be  linearly  dependent  upon  temperature 
for 

4500  'temperatures  greater  than  about  73.5  degrees. 

4510  ' 

4520  ' 

4530  TF=  TEMP-73.5 

4540  IF  TF>0  THEN  READING (3) =  READING (3) +. 0565*TF 

4550  'TF  is  used  as  a  variable  name  for  temperature  difference 

4560  ' 

4570  ' 

4580  ' 

4590  EAA=READING (2) -OAA  'The  offsets  are  applied  to  the 
channels 

4600  EBA«READING(3)-OBA 
4610  EAN»READING(4)-OAN 
4620  EBN»READING{5)-OBN 
4630  ' 

4640  ' 

5130  ' 

5140  ' 

5150  ' 

5160  ' 

5170  FOR  I  =  1  TO  10 
5180  ' 

5190  ' 

5500  ' 
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5510  ' 

5520  ' 

5530  t  *****************  Balance  Interaction  Equations 

ic'kieitititit'k'kic'kiticldtit'k 

5540  ' 

5550  'These  balance  interaction  equations  were  developed  by 
5560  'John  Kersh.  You  must  ensure  that  they  still  apply  via 
a 

5570  'calibration  procedure.  The  calibration  procedure 
involves 

5580  'hanging  known  weights  off  the  calibration  rig  at  known 
turn 

5590  'table  angles  of  attack.  The  actual  normal  and  axial 
forces 

5600  'are  then  found  using  Normal®  Weight *cos (AOA)  and  Axial® 
5610  ' Weight *sin (AOA) . 

5620  ' 

5630  ' 

5640  ' 

5650  ' 

5660  MAN®27 .26281 *EAN  -  .487970*EAA  'Watch  sign  convention 
here 

5670  MAA®-.0493204*EAN  +  32 .538138*EAA 
5680  MBN®11.303456*EBN  -  .7095687*EBA 
5690  MBA®-. 1350 659*EBN  +  17 .258483*EBA 

5700  BB=2.2083  'Separation  of  upper  and  lower  strain  bridges, 
5710  '  'B  and  A,  in  feet. 

5720  'Large  parts  of  this  program  were  untouched! 

5730  'A  marker  for  the  iterations 
5740  CYCLE  ®  CYCLE  +  1 
5750  NEXT  I 
5760  ' 

5770  ****************  force  EQUATIONS  ******************* 
5780  ' 

5790  ' 

5800  'Normal,  Axial,  Lift,  and.  Drag  forces  are  calculated 
here. 

5810  ' 

5820  ' 

5830  NORMAL  ®  (MAN-MBN) /BB 
5840  ' 

5650  AXIAL  ®  (MAA-MBA)/BB 

5  8  6  0  DRAG®  A X  I  A L * S I N  (  .  0  1  7  4  5  3  2  9  *  A 0 A ) 
NORMAL*COS ( . 01745329*AOA) 

5870  LIFT®  AXIAL*COS ( . 01745329*AOA)  + 

NORMAL*SIN ( .01745329*AOA) 

5880  ' 

5890  ' 

5900  TRIAL  ®  TRIAL  +  1 

5910  INPUT  "THIS  IS  NOT  A  TARE  READING.  HIT  <CR>";AN$ 

5920  IF  AN$  <>  "Y"  GOTO  6190 


5930  COLOR  0,10,10:CLS 

5935  IF  A$="Y'’  THEN  TRIAL=TRIAL-'l  :GOTO  5950 
5940  TRIAL  =  0 
5950  TRY  »  TRY  +  1 
6000  TARE(l)  *  TRY 
6010  TARE  (2)  =  AOA 
6020  TARE (3)  =  NORMAL 
6030  TARE (4)  «  AXIAL 
6040  TARE (5)  =  LIFT 
6050  TARE (6)  =  0 
6060  TARE (7)  =  DRAG 
6070  TARE (8)  =  0 
6080  ' 

6090  '  PRINT  THE  TARING  DATA 
6100  PRINT" 

ffl 


6110  PRINT"*  *  *  * 

***** 

*  *  * 

TARE  CALCULATIONS 

it  ik  ik  it 

6120  PRINT 

"  TRIAL 

AOA 

NORMAL 

SIDE 

AXIAL 

PITCH  ROLL 

6130  PRINT 

YAW" 

"  # 

DEG 

POUNDS 

POUNDS 

POUNDS 

FT-LBS  FT-LBS 

6140  PRINT 

FT-LBS" 

n  *****  ***** 

****** 

****** 

****** 

******  ****** 

6150  ' 

6160  PRINT  USING" 

******** 

##  +## 

.«« 

«««.«« 

«««.«« 

«««.«« 

«  «  «  .  # 

# 

«  «  « 

# 

« 

### . ##";TARE (1) , TARE (2) , TARE (3) ,TARE (4) , TARE (5) , TARE (6) , TARE 

(7)  ,TARE(8) 

6170  BEEP: INPUT  "ENTER  <CR>  TO  CONTINUE"; INPT$ : IF  A$="Y"  THEN 
GOTO  6350 
6180  GOTO  6460 
6190  ' 

6240  FORCE (TRIAL, 1)  =  TRIAL 

6250  FORCE (TRIAL, 2)  =  AOA 

6260  FORCE (TRIAL, 3)  »  CNDAOA 

6270  FORCE (TRIAL, 4)  *  NORMAL 

6280  FORCE (TRIAL, 5)  «  AXIAL 

6290  FORCE (TRIAL, 6)  ■  LIFT 

6300  FORCE (TRIAL, 7)  »  DRAG 

6310  FORCE (TRIAL, 8)  «  TEMP 

6320  FORCE (TRIAL, 9)  »  TEMP 

6330  'print  the  values  and  store  in  file 

6340  ' 

6350  PRINT" 

tf 

6360  PRINT"*  ***********  FORCE  CALCULATIONS  *  *  *  * 
********* 
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6370  PRINT  -  TRIAL  AOA  CANARD  NORMAL  AXIAL  LIFT 
DRAG  TEMP" 

6371  'YAW" 

6380  PRINT 

6381  PRINT  "  #  DEG  DEG  LBS  LBS  LBS 

LBS  F" 

6390  'PRINT 

6391  '  *****  *****  ******  ******  ******  ****** 

******  ******" 

6400  'a  loop  to  list  all  values  so  far 
6410  ' 

6420  FOR  J  =  1  TO  TRIAL 

6430  PRINT  USING"  ##  +###.#  -<-###.##  +###.##  +###.## 
+  ###.##  +###.## 
+###.##"; FORCE ( J, 1 ) , FORCE (3,2), FORCE ( J, 3 ) , FORCE ( J, 4 ) , FORCE ( J 
, 5) , FORCE ( J, 6) , FORCE ( J, 7) , FORCE (J, 8) 

6440  NEXT  J 

6450  BEEP: INPUT  "ENTER  <CR>  TO  CONTINUE"; INPT$ 

6460  ' 

6470  'Writes  the  data  to  the  output  data  files 

6480  COLOR  14,1,1:CLS 

6490  OPEN  DATAFILES  FOR  OUTPUT  AS  #1 

6500  OPEN  DISKFILES  FOR  OUTPUT  AS  #2 

6510  WRITE#1, 

TARE (1) , TARE (2) , TARE (3) , TARE (4) , TARE (5) , TARE (6) , TARE (7) , TARE 
(8) 

6520  WRITE#2, 

TARE  (1) ,  TARE  (2) ,  TARE  (3) ,  TARE  (4)  ,  TARE  (5)  ,  TARE  (6)  ,  TARE  (7)  ,  TARE 
(8) 

6530  FOR  X  *  1  TO  140 

6550  WRITE  #1. 

FORCE (X, 1 ) , FORCE (X, 2) , FORCE (X, 3) , FORCE (X, 4 ) , FORCE (X, 5) , FORCE 
(X,  6) , FORCE (X, 7) , FORCE (X, 8) , FORCE (X, 9) 

6560  WRITE  #2, 

FORCE (X,l) , FORCE (X,2) , FORCE (X, 3) , FORCE (X,4) , FORCE (X, 5) , FORCE 
(X,  6)  , FORCE (X, 7) , FORCE (X, 8) , FORCE (X, 9) 

6570  NEXT  X 
6580  CLOSE  #1 

6590  CLOSE  #2  'To  view  the  corrected  voltage  output 

6600  LOCATE  9,10:PRINT  "THE  CORRECTED  CHN(2)  READING  WAS";EAA 
6610  LOCATE  10,10:PRINT  "THE  CORRECTED  CHN(3)  READING  WAS";EBA 
6620  LOCATE  11,10:PRINT  "THE  CORRECTED  CHN(4)  READING  WAS";EAN 
6630  LOCATE  12,10:PRINT  "THE  CORRECTED  CHN(3)  RE?d)ING  WAS";EBN 
6640  LOCATE  14, 10: INPUT  "DO  YOU  WANT  ANOTHER  SCAN  (Y  OR 
N)";ANSW$ 

6645  A$*"N" 

6650  IF  ANSW$  <>"N"  THEN  GOTO  3700 

6655  CLS: LOCATE  12, 15 .-INPUT  "DO  YOU  REALLY  WANT  TO  QUIT";AW$ 
6660  IF  AW$<>"Y"  THEN  GOTO  3700 
6665  60SUB  7060 
6670  END 
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6680  ' 

6690  'This  subroutine  averages  the  balance  voltage  readings 
6700  'by  computing  the  mean  and  standard  deviation. 

6710  'Any  readings  less  or  greater  than  one  standard  deviation 
6720  'are  thrown  out  and  a  new  mean  is  computed 
6730  ' 

6740  FOR  CHANNEL  =  2  TO  5 
6750  N=10:FLAG=0 
6760  SSDEV=0 

6770  'Mean  of  balance  voltage  readings 

6780  SREAD  »  0 

6790  FOR  CNT  *  1  TO  10 

6800  SREAD  *  SREAD  +  TREAD (CHANNEL, CNT) 

6810  NEXT  CNT 

6820  MEAN  =  SREAD/N 

6830  READING (CHANNEL)  «  MEAN 

6840  IF  {FLAG=1)  THEN  GOTO  7010 

6850  'Standard  deviation  routine 

6860  FOR  CNT  *  1  TO  10 

6870  DIF  =  TREAD  (CHANNEL,  CNT)  -  MEAN 

6880  SDEV  *  DIF  *  DIF 

6890  SSDEV  *  SSDEV  +  SDEV 

6900  NEXT  CNT 

6910  DEV  =  SQR(SSDEV/N) 

6920  HI  =  MEAN  +  DEV 
6930  LO  =  MEAN  -  DEV 
6940  FOR  CNT  *  1  TO  10 
6950  ARG  =  TREAD (CHANNEL, CNT) 

6960  IF  (ARG  <  HI)  AND  (ARG  >  LO)  THEN  GOTO  6990 
6970  TREAD  (CHANNEL,  CNT)  »=  0 
6980  N  »  N  -  1:FLAG=1 
6990  NEXT  CNT 

7000  DROP (CHANNEL) =10-N: GOTO  6780 
7010  NEXT  CHANNEL 

7020  PRINT  "READINGS  DROP  ";DROP(2);"  ";DROP(3);" 

";DROP(4);"  ";DROP(5);"  ";DROP(6);" 

"/DROP  (7) 

7030  RETURN 
7040  END 
7050  ' 

7060  'THIS  SUBROUTINE  ENTER  THE  EXPERIMENT  CONDITIONS. 

7070  COLOR  14,5,5:CLS 

7080  LOCATE  12, 10  .'INPUT  "ENTER  EXPERIMENT  DATE 
(YYMMDD)  "  ; YMD : LAB  ( 1 )  «YMD 

7090  LOCATE  13, 10: INPUT  "ENTER  THE  EXPERIMENT  GRID 
NO.";G:LAB(2)»G 

7100  LOCATE  14,10:INPUT  "WAS  THE  CANARD  IN  PLACE?  l»NO, 
2»YES";B:LAB(3)-B 

7110  LOCATE  15, 10: INPUT  "ENTER  THE  NEXT  STARTING 

AOA";S:LAB(4)«S 

7115  TEMPSUM=0 


7120  FOR  X*1  TO  TRIAL 

7125  TEMPSUM=TEMPSUM+FORCE(X, 9) 

7130  NEXT  X 

7135  TAV*TEMPSUM/TRIAL : LAB ( 5 ) =TAV 

7140  LOCATE  17, 10: INPUT  "ENTER  EXPR  TUNNEL  DYN  PRESSURE 
(cmH20) ";DP:LAB(6)=DP 

7150  LOCATE  18, 10: INPUT  "ENTER  EXPR  PRESSURE  (in. 

Hg) ";PRE:LAB(7)*PRE*70.739 
7160  ' 

7170  OPEN  DATAFILES  FOR  OUTPUT  AS  #1 
7180  OPEN  DISKFILES  FOR  OUTPUT  AS  #2 

7190  WRITE#1, 

TARECl)  ,  TARE  (2)  ,  TARE  (3) ,  TARE  (4) ,  TARE  (5)  ,  TARE  (6) ,  TARE  (7) ,  TARE 

7200  WRITE#2, 

TARE ( 1) , TARE (2) , TARE (3) , TARE (4) , TARE (5) , TARE (6) , TARE (7) , TARE 
(8) 

7210  FOR  X  =  1  TO  140 

7220  WRITE 

#1 , FORCE (X, 1 ) , FORCE (X, 2 ) , FORCE (X, 3) , FORCE (X, 4 ) , FORCE (X, 5) , FO 
RCE (X, 6) , FORCE (X, 7 ) , FORCE (X, 8 ) , FORCE (X, 9 ) 

7230  WRITE  #2, 

FORCE (X, 1) , FORCE (X, 2) , FORCE (X, 3) , FORCE (X, 4) , FORCE (X, 5) , FORCE 
(X, 6) , FORCE (X, 7) , FORCE (X, 8) , FORCE (X, 9) 

7240  NEXT  X 

7250  WRITE#1, 

LAB  ( 1 ) ,  LAB  (2 ) ,  LAB  (3) ,  LAB  (4 ) ,  LAB  (5) ,  LAB  ( 6) ,  LAB  (7 ) 
7260  WRITE#2, 

LAB  (1)  ,  LAB  (2)  ,  LAB  (3) ,  LAB  (4) ,  LAB  (5) ,  LAB  (6) ,  LAB  (7) 

7270  CLOSE  #1 
7280  CLOSE  #2 
7290  ' 

7300  ' 

7310  CLS: LOCATE  12,10 

7320  PRINT  "GAME  OVER!  MAN! " :BEEP : BEEP 
7330  RETURN 
7340  END 


M 
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APPENDIX  D.  DATA  MANIPULATION/CORRZCTION  PROGRAM 


1000  •  PROGRAM  BY  M.P.  RABANG  AND  LT  YUAN,  C.C  TO  READ  FORCE 
AND 

1005  '  MOMENT  VALUES  FROM  A  DATA  AQUISITION  PGM  TITLED 

READ. BAS 

1010  •  THIS  PROGRAM  WAS  MODIFIED  LT  J.M.  KERSH  TO  MANIPULATE 
THE 

1015  '  THE  PGM  ENTITLED  JMK.BAS.  JMK.BAS  IS  THE  DATA 
AQUISITION 

1020  •  PROGRAM  FOR  THE  CLOSE-COUPLED  CANARD /WING  MODEL. 

1030  DIM  TARE[8],  FORCE [140, 9) /  COEFI140,9],  LAB[11] 

1040  COLOR  14,1,1 
1050  CLS 
1060  KEY  OFF 

1070  /**************Last  Modified  on  18  OCT  90 

**4t*****ik*********** 

1080  LOCATE  11,7 

1090  INPUT"ENTER  THE  NAME  OF  THE  INPUT  FILE";D$ 

1100  F$=D$+".DAT" 

1110  INFILE$»"C;\KERSH\"+F$ 

1120  CF$»"CF"+F$ 

1130  OUTFILE$=''C;\KERSH\"+CF$ 

1140  DISKFILE$=''A:"+CF$ 

1150  '  Force  values  are  read  from  the  input  data  file 
1160  OPEN  INFILES  FOR  INPUT  AS  #1 

1170INPUT#1, 
TARE (1) , TARE (2) , TARE (3) , TARE (4) , TARE (5) , TARE (6) , TARE (7) , TARE 
(8) 

1180  FOR  X  =  1  TO  140 

1190INPUT#1, 
FORCE (X, 1) , FORCE (X, 2) , FORCE (X, 3) ,FORCE(X,4) ,FORCE(X,5) , FORCE 
(X,  6) ,  FORCE  (X,  7) ,  FORCE  (X,  8) ,  FORCE  (X,  9) 

1200  NEXT  X 

12  10INPUT#1, 
LAB  (1) ,  LAB  (2) ,  LAB  (3) ,  LAB  (4) ,  LAB  (5) ,  LAB  (6)  ,  LAB  (7) 

1220  CLOSE  «1 
1225  ' 

1230  YMD=LAB(1)  'YEAR  MONTH  DATE 

1240  SCR-LAB(2)  'EXPERIMENT  GRID  NUMBER 

1250  BODY-LAB (3)  'WAS  THE  CANARD  IN  PLACE 

1260  STR-LAB(4)  'ANGLE  OF  ATTACK 

1270  TAV-LAB(5)  'AVERAGE  TEMPERATURE 

1280  DP-LAB (6)  'DYNAMIC  PRESSURE 

1290  PRE-LAB (7)  'ATMOSPHERIC  PRESSURE 

1300  ' 

1330  '  OPEN  THE  DATA  FILE  SO  EACH  SCAN  IS  RECORDED 
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1340  TRANSFILE$="C: \KERSH\TRANS.DAT" 

1370  ' 

1400  CLS: LOCATE  12,5 

1410  INPUT"  INPUT  FILE  HAS  BEEN  LOADED,  ENTER  <CR>  TO 
CONTINUE" ;INPT$ 

1420  ' 

1430  ^ **************  BEGIN  COEFFICIENTS  CALCULATION 

************ 

1440  IF  BODy=l  THEN  A®. 681  '  Area  Of  wing  to  centerline  of 

fuselage 

1450  MU®3.719E-07  • 

1460  RH0=PRE/ (1545* (459. 7+TAV) ) 

1470  Q®-. 026749+1. 1149*DP 

1480  IF  BODY®2  THEN  A®. 815  'Area  of  wing  and  exposed  area  of 
canard 

1490  Q*Q*2.047 
1500  VEL=SQR(2*Q/RHO) 

1510  RED®  (RHO*VEL*  (9.52/12)) /MU  'Ref  Area  is  MAC  of  Wing  to 

centerline 

1520  FLAG®0 

1530  FOR  X  =  1  TO  140 

1540  IF  FORCE (X,1)=0  THEN  GOTO  1730 

1550  FLAG®FLAG+1 

1560  COEF(X,  DEFORCE  (X,  1) 

1570  COEF(X,2)=FORCE(X,2)-2.2  'Corrects  for  balance 

misalignment 

1580  COEF(X,3)«FORCE(X,3) 

1590  COEF(X,4)®.955*FORCE(X,4)  'Corrects  for  calerror  in 
normal  force 

1600  COEF(X,5)=1.037*FORCE(X,5)  'Corrects  for  calerror  in 
axial  force 

16  10 

COEF (X, 6) ®COEF (X,  5) *COS ( . 017453*COEF (X, 2) ) +COEF (X, 4) *SIN ( . 01 
7453*COEF  (X,2) ) 

16  2  0 

COEF (X, 7) =COEF (X,  5) *SIN ( . 017453*COEF (X, 2) ) -COEF (X, 4) *COS  ( .  01 
7453*COEF  (X,2)) 

1630  f  ************  Blockage  Correction  Routine 

********************** 

1640  IF  COEF(X,2)  =  0  THEN  EPS®. 00269  'Area  of  4.5x3  inch 
block 

1650  EPBW=.0298*SIN(.017453*COEF(X,2) ) 

1660  EPCND®.00383*SIN  ( (COEF (X, 2) +COEF (X, 3) ) * . 017453) 

1670  IF  BODY«l  THEN  EPS=ABS (EPBW) 

1680  IF  BODY«2  THEN  EPS»EPBW+EPCND 

1700  COEF  (X, 9) -COEF (X, 7) / (A*Q* (1+ (2*EPS) ) ) 

1710  COEF(X,8)»COEF(X,  6)/(A*Q*(l+(2*EPS))) 

1720  'Original  pgm  CORR>BAS  greatly  changed  to  account  for 
turntable 

1730  'misaligiunent .  Turntable  was  +2.2  degrees  off. 

1740  NEXT  X 
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1750  LAB(8)*RHO 
1760  LAB(9)=Q 
1770  IiAB(10)=VEL 
1780  LAB(11)=RED 

1790  OPEN  OUTFILE$  FOR  OUTPUT  AS  #1 
1800  OPEN  DISKFILE$  FOR  OUTPUT  AS  #2 
1810  OPEN  TRANSFILE$  FOR  APPEND  AS  #3 

1820  WRITE  #1, 

LAB  (1) ,  LAB  (2)  ,  LAB  (3)  ,  LAB  (4)  ,  LAB  (5)  ,  LAB  (6)  ,  LAB  (7) ,  LAB  (8)  ,  LAB  ( 
9)  ,  LAB  (10)  ,  LAB  (11) 

1830  WRITE  #2, 

,  LAB(l)  ,LAB(2)  ,LAB(3)  ,LAB(4)  ,  LAB  (5)  ,  LAB  (6)  ,  LAB  (7)  ,  LAB  (8)  ,  LAB  ( 

9)  ,  LAB  (10) ,  LAB  (11) 

1840  WRITE  #3, 

LAB(l)  ,LAB(2)  ,LAB(3)  ,LAB(4)  ,LAB(5)  ,LAB(6)  ,LAB(7)  ,LAB(8)  ,LAB( 
9) ,  LAB  (10)  ,  LAB  (11) 

1850  FOR  X=1  TO  FLAG 

1860  WRITE#!  ', 

COEF  (X,  1)  ,  COEF  (X,  2) ,  COEF  (X,  3)  ,  COEF  (X,  4) ,  COEF  (X,  9) ,  COEF  (X,  5)  , 
COEF (X, 6) , COEF (X, 7 ) , COEF (X, 8 ) 

1870  WRITE#2, 

COEF(X,l)  ,COEF(X,2)  ,COEF(X,3)  ,COEF(X,4)  ,COEF  (X,  9)  ,  COEF  (X,  5)  , 
COEF (X, 6) , COEF (X, 7) , COEF (X, 8) 

1880  WRITE#3, 

COEF(X,  1)  ,COEF(X,2)  ,COEF(X,3)  ,COEF(X,4)  ,COEF  (X,  9)  ,  COEF  (X,  5)  , 
COEF (X, 6) , COEF  (X, 7) , COEF (X, 8) 

1890  NEXT  X 
1900  CLOSE  #1 
2000  CLOSE  #2 
2010  CLOSE  #3 
2020  ' 

2060  •  DISPLAY  ROUTINE 

2070  CLS; BEEP: LOCATE  10, 5: INPUT  "DO  YOU  WANT  TO  VIEW  THE 
OUTPUT"; PANS $ 

2080  IF  PANS$<>"Y"  THEN  GOTO  2290 
2090  COLOR  0,10,10 
2100  CLS 

2110  PRINT"FILE  NAME:  ";CF$: PRINT"  " 


^  2120  PRINT "DATE  (YYMMDD)  ";YMD 

2125  PRINT" SCREEN  NO.  ";SCR 

2130  PRINT"WAS  THE  CANARD  IN  PLACE?  l»NO,  2«YES  ";BODY 
‘  2135  PRINT"THE  LAST  AOA  FOR  THE  RUN  WAS  ";STR 

2140  PRINT"STATIC  PRESSURE  (LB/FT^2)  ";PRE 

2145  PRINT"AVERAGE  TEMPERATURE  (F)  ";TAV 

2150  PRINT "WIND  TUNNNEL  VELOCITY  (FT /SEC)  ";VEL 

2155  PRINT"WIND  TUNNNEL  DYNAMIC  PRESSURE  (cinH20)  ";DP 

2160  PRINT"AIR  DENSITY  (LBm/FT''3)  ";RHO 

2170  PRINT"REYNOLDS  NUMBER  ";RED 

2180  PRINT"ACTUAL  DYNAMIC  PRESSURE  (LB/FT^2)  ";Q 

2185  BEEP -.INPUT  "ENTER  <CR>  TO  CONTINUE";  INPT$ 

2190  PRINT"  " 
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2200  PRINT"*  *************  FORCE  COEFFICIENTS  * 


********* 

2210  PRINT"  " 

2220  PRINT  "TRIAL 

*  *  *  * 

AOA  CNDAOA 

NORMAL 

AXIAL 

LIFT 

DRAG  TEMP  " 

2221  ' 

2225  ' 

2230  PRINT  "***** 

***  ****** 

****** 

***** 

*  *  *  * 

****  ****n 

2240  FOR  X=1  TO  FLAG 

2250  PRINT  USING  "  ###  +###.# 

+##.# 

+##.### 

•>■««.««« 

+  ««.«« 

#  +  # 

«  .  « 

«  « 

##.##";COEF  (X,l),. 

COEF (X,2), COEF (X, 

3),  COEF  (X 

,4)  ,COEF(X 

i,  5),COE 

F  (X,  6) ,  COEF  (X,  7 ) ,  COEF  (X,  8 ) 

2255  IF  X=20  OR  X=40  OR  X=60  OR  X*80  THEN  INPUT  "ENTER  <CR>  TO 
CONTINUE" ;INPT$ 

2257  ' 

2258  ' 

2260  ' 

2270  NEXT  X 

2280  BEEP  .-INPUT  "ENTER  <CR>  TO  CONTINUE";  INPT$ 

2290  CLS: LOCATE  1 0 ,  5 : BEEP :  INPUT  "DO  YOU  WANT  A  HARDCOPY"  ;ANS$ 
2300  IF  ANS$<>"Y"  THEN  GOTO  2640 

2310  LPRINT  CHR$ (27) ;"&16E"  'Output  changed  to  thesis 

format 
2312 
2314 
2320 
2345 
2350 
2355 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 

****** 

2460  LPRINT"  " 

2470  LPRINT"TRIAL 
DRAG  TEMP" 

2480  LPRINT"  # 

DEG  F" 

2490  LPRINT"***** 

****  *****" 

2500  FOR  J  *  1  TO  FLAG 

2520  LPRINT  USING"  ###  +###.#  +##.#  +##.###  +##.###  +##.### 


AOA  CNDAOA  NORMAL  AXIAL  LIFT 
DEG  DEG  LBS  LBS  LBS  LBS 


*** 


*** 


**** 


*  *  *  * 


**** 


LPRINT  CHR$ (27) ; "&149F"  'These  are  toggles 
LPRINT  CHR$ (27) ; "ialOL" 

LPRINT"FORCE  DATA  FILENAME: ";F$ 

LPRINT"COEFFICIENT  DATA  FILENAME: ";CF$ : LPRINT"  " 
LPRINT"DATE  (YYMMDD)  ";YMD 

LPRINT"WAS  THE  CANARD  IN  PLACE?  l=NO, 2*YES  ";BODY 

LPRINT "ATMOSPHERIC  PRESSURE  (Inches  Hg)  ";PRE/70.739 
LPRINT "AVERAGE  TEMPERATURE  (F)  ";TAV 

LPRINT"WIND  TUNNEL  VELOCITY  (FT/SEC)  ";VEL 

LPRINT"WIND  TUNNEL  DYNAMIC  PRESSURE  (cmH20) ";DP 
LPRINT"AIR  DENSITY  (LBM/FT''3)  ";RHO 

LPRINT"REYNOLDS  NUMBER  ";RED 

LPRINT"ACTUAL  DYNAMIC  PRESSURE  (LB/FT"2)  ";Q 

LPRINT"  " 

LPRINT"  " 

LPRINT"*  **********  force  READINGS  *  *  *  *  * 


116 


+  ##.### 

##.##"; FORCE (J,l), FORCE (J,2), FORCE ( J, 3) , FORCE ( J, 4 ) , FORCE ( J, 5 
), FORCE  (J,  6)  ,FORCE(J,7)  ,FORCE(J,  8) 

2530  NEXT  J 
2540  LPRINT"  " 

2550  LPRINT"  " 

2560  LPRINT"*  *  *  *  *  CORRECTED  AOA,  CL,  CD,  LIFT,  AND  DRAG 
***** 

2570  LPRINT"  " 

2580  LPRINT"TRIAL  AOA  CNDAOA  C  LIFT  C  DRAG  LIFT 

DRAG  " 

2590  LPRINT"*****  ***  ***  ****  ****  **** 

****  1 

2600  FOR  X  *  1  TO  FLAG 

2610  LPRINT  USING"  ###  +###.#  +##.#  ->-««. «««#  +##.#### 

+  ##.####  +§§.*§** 

";COEF (X, 1) , COEF (X, 2) , COEF (X, 3) , COEF (X, 8) , COEF (X, 9) , COEF (X, 6 
)  ,COEF(X,7) 

2630  NEXT  X 

2640  COLOR  14, 1, 1 :CLS : LOCATE  12,7:BEEP 
2650  INPUT"DO  YOU  WANT  ANOTHER  RUN";AANS$ 

2660  IF  AANS$<>"N"  THEN  GOTO  1050 

2670  CLS; LOCATE  12, 10: PRINT  "GAME  OVER!  MAN !": BEEP : BEEP 
2680  END 


APPENDIX  Z.  DATA  SETS 


Table  6.  Baseline  Configuration  and  Comparison  to 
Canard/Wing/Body;  Part  A 


I 


■11.1700 


-7.9040 


-4.9937 


-2.0504 


.9644 


4.0655 


7.4842 


11.7310 


zssm\ 


17.9998 


man\ 


22.7860 


25.0035 


26.6030 


27.7366 


28.3916 


29.4518 


28.6064 


28.5706 


28.6415 


28.6372 


30.6738 


31.1688 


33.6451 


34.6974 


38.7898 


40.9079 


2.0010 


1.5016 


1.2056 


0.9030 


0.8945 


1.1027 


1.5756 


2.2435 


3.1626 


4.1417 


5.4213 


6.9427 


8.5957 


9.7539 


10.4862 


43.0890 


43.4319 


43.0571 


41.2888 


39.7179 


11.8033 


12.4834 


13.2223 


13.8457 


15.2944 


16.1211 


18.1506 


19.2472 


20.8344 


23.6592 


26.6938 


29.3636 


32.1095 


34.5606 


36.7935 


37.5997 


36.3503 


■KIKIEQ  ■KIEIB] 


a 


0.4198 


.2977 


0.1885 


0.0775 


.0365 


0.1538 


.2826 


.4420 


.5555 


0.6754 


.7618 


0.6515 


.9326 


0.9903 


1.0548 


1.0932 


1.0608 


1.0584 


1.0601 


1.0589 


1.1332 


1.1504 


1.2407 


1.2784 


1.3458 


1.4254 


1.5007 


1.5483 


1.5756 


1.5857 


1.5696 


1.5029 


1.4437 


lESMin 


CD 


0.0752  23AUGA.DAT 


.0566 


0.0454 


0.0388 


0.0342 


0.03  3  8  23AUGA.DAT. 


0.0416 


.0594 


.0844 


0.1187 


.155  1  23AUGB.DAT 


0.2026  B 


0.2589  B 


20SEPA.DAT 


0.3896 


0.4249 


0.4377 


.4625 


0.4894 


.5120 


.5650 


.5950 


.6693 


.7091 


.7670  23AUGC.DAT 


0.6694  C 


.9793IC 


C 


1.1741  23AUGD.DAT 


1.2625 


37.5366 


40.6866 


1.3608 


1.3413 


1.3687 


1.3940 


1.4308 


1.4750 


Table  7.  Baseline  Configuration  and  Comparison  to 
Canard/Wing/Body;  Part  B 
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Table  8.  Baseline  Configuration  and  Comparison  to 
Canard/Wing/Body;  Part  C 
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Table  9.  Wing/Body  at  a  10  Degree  AOA 


1  CANARDAOA 

maoQsi 

a 

CD 

DATA  RUN 

15.9395 

3.6655 

0.5007 

0.1152 

20SEPC.DAT 

"IWTTll 

3.2091 

0.5200 

0.1007 

C 

-15 

3.0992 

0.5675 

0.0972 

C 

-10 

19.2692 

3.1160 

0.6041 

0.0977 

C 

-  5 

20.1266 

3.1740 

0.6306 

0.0944 

C 

0 

21.2552 

3.4418 

0.6655 

0.1078 

23SEPG.DAT 

5 

21.1985 

3.8438 

0.6633 

0.1203 

G 

7 

21.8989 

4.2309 

0.6851 

0.1324 

G 

1  0 

21.6258 

4.5604 

0.6763 

0.1426 

G 

12 

22.2098 

4.9784 

0.6944 

0.1556 

G 

1  5 

5.5574 

0.6962 

0.1737 

G 

1  7 

22.3500 

6.0144 

0.6983 

0.1879 

G 

20 

22.2796 

6.5067 

0.6959 

0.2032 

G 

25 

21.2853 

0.6645 

0.2227 

G 

Table  10.  Wing/Body  at  a  22  Degree  AOA 


CANARD  AOA 

DRAG  (LBS) 

a 

CD 

DATA  RUN 

-25 

31.3361 

11.8939 

0.9713 

0.3687 

22SEPA.DAT 

-20 

32.5169 

11.8803 

1.0073 

0.3680 

A 

-15 

12.1208 

1.0702 

0.3752 

A 

-10 

36.4250 

12.7764 

1.1269 

0.3953 

A 

-  7 

37.7855 

13.3625 

1.1685 

0.4132 

A 

-  5 

39.3184 

13.8677 

0.4288 

23SEPA.DAT 

-  2 

41.7156 

14.8393 

1.2893 

0.4586 

A 

0 

44.2310 

15.9831 

1.3667 

0.4939 

A 

2 

44.9294 

16.2654 

1.3879 

0.5025 

A 

5 

45.7376 

16.9976 

1.4124 

0.5249 

A 

7 

46.0463 

17.7145 

1.4216 

0.5469 

A 

1  0 

■Kmoi 

18.1777 

1.4109 

0.5610 

A 

12 

43.6964 

17.8816 

1.3483 

0.5518 

A 

15 

43.4425 

18.3095 

1.3400 

0.5648 

A 

20 

43.2672 

19.0409 

1.3340 

0.5870 

A 

25 

42.8524 

19.6500 

1.3206 

0.6055 

A 
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Table  11.  Wing/Body  at  a  34  Degree  AOA 


DRAG  (LBS) 

a 

CD 

DATA  RUN 

-25 

44.5152 

25.3760 

1.3630 

0.7770 

17SEPC.DAT 

-20 

■Easi 

— iwnTi 

1.4076 

0.8027 

C 

-15 

48.1965 

27.4527 

1.4739 

0.8395 

C 

-10 

51.5752 

30.0208 

1.5763 

0.9175 

-  7 

32.1670 

1.6420 

0.9828 

A 

-  5 

52.9515 

32.1368 

1.6174 

0.9816 

A 

-  2 

52.6874 

32.7322 

1.6088 

0.9995 

A 

0 

32.7426 

1.5913 

0.9996 

A 

2 

52.4038 

33.5410 

1.5995 

1.0237 

A 

5 

53.2707 

34.6841 

1.6254 

1.0583 

A 

7 

52.4110 

34.4708 

1.5989 

1.0516 

25SEPB.DAT 

to 

51.9213 

34.8199 

1.5835 

1 .0619 

B 

1  S 

51.7202 

35.5922 

1.5766 

1.0850 

B 

20 

51.2194 

36.1678 

1.5607 

1.1021 

B 

25 

TFnn 

36.5055 

1.5387 

1.1120 

B 

Table  12.  Wing/Body  at  a  40  Degree  AOA 


CANARD  AOA 

DRAG  (LBS) 

CL 

CD 

DATA  RUN 

•25 

52.1035 

36.0908 

1.5865 

1.0989 

18SEPD.DAT 

-20 

54.1898 

37.9837 

1.6490 

1.1559 

18SEPB.DAT 

-17 

54.9179 

38.7684 

1.6706 

1.1793 

24SEPB.DAT 

-15 

55.8909 

39.9075 

1.6998 

1.2137 

B 

-10 

55.8655 

40.0806 

■BSil!] 

■nfTFFi 

B 

-  7 

55.4968 

41.0806 

1.6863 

sbheilli 

B 

-  5 

55.6604 

41.5648 

1.6909 

1.2627 

B 

-  2 

55.3475 

42.2035 

1.6609 

1.2817 

B 

0 

55.3350 

42.7420 

1.6801 

■HiSEQ 

B 

5 

54.7972 

43.6663 

1.6630 

1.3252 

18SEPA.DAT 

10 

■TTirm 

44.1921 

1.6382 

1.3406 

A 

15 

52.3813 

43.9563 

1.5884 

1.3330 

18SEPD.DAT 

20 

51.4276 

44.2213 

1.5589 

1.3330 

D 

25 

49.9234 

43.9899 

1.5129 

1.3405 

0 
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) 


Table  13.  Wing/Body  at  a  48  Degree  AOA 


1  CANAROAOA 

a 

CD 

DATA  RUN 

63.6236 

46.7047 

1.6219 

1.4731 

25SEPADAT 

54.0065 

49.7988 

1.6329 

1.5057 

A 

54.1415 

49.8229 

1.6366 

1.5060 

A 

-17 

wtmEm 

50.7786 

1.6491 

1.5344 

A 

•  IS 

54.1012 

51.0222 

1.6345 

1.5414 

A 

-10 

53.5873 

52.1664 

1.6181 

1.5752 

A 

-  5 

53.3432 

53.2749 

1.6099 

1.6079 

19SEPA.DAT 

0 

52.6276 

53.7827 

1.5876 

1.6225 

19SEPB.DAT 

5 

51.1209 

53.2744 

1.5416 

1.6065 

B 

1  0 

49.7827 

53.1249 

1.6014 

B 

15 

48.1900 

52.5577 

B 

20 

46.5715 

51.8891 

1.4030 

1.5632 

B 

25 

45.7369 

52.2085 

1.3776 

B 

t 
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